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ABSTRACT 

' Beta-decay and muon-capture experiments in nuclei are reviewed. The conserved vec- 

tor current hypothesis is confirmed through the observed constancy of the vector coupling 
5-H , constant determined from the superallowed Fermi transitions and from the measurement 

of the weak-magnetism term in mirror Gamow- Teller transitions. The axial-vector and 
pseudoscalar coupling constants in the nucleon are determined from neutron decay and 
muon capture on the proton respectively. In finite nuclei, evidence for these coupling con- 
stants being reduced relative to their free-nucleon values is discussed. Meson-exchange 
currents are shown to be an important correction to the time-like part of the axial cur- 
rent as evident in first-forbidden beta decays. Tests of the Standard Model are discussed, 
as well as extensions beyond it involving right-hand currents and scalar interactions. 
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1.1. Standard Model for Quarks and Leptons 

The standard electroweak model is based on the gauge groupS SU{2) x ?7(1), 
with four vector gauge bosons VF^, a = 1, 2, 3, and for the SU{2) and U{1) sectors 
^ ! respectively, and the corresponding gauge coupling constants g and g' . The left-handed 

\ fermion fields, = (vi^i^) and {ui,d'^), comprise neutrinos, i^j, leptons, i^, up and 

down quarks, Ui and d'^ for the i^^ fermion family, and transform as doublets under 



X 



SU{2), where d[ = J2j ^ijdj ^ is the Cabibbo-Kobayashi-Maskawa (CKM) matrix, 
which we discuss shortly. The right-handed fields are SU{2) singlets. The interaction 



! part of the Lagrangian for quarks and leptons is eII 



IP 

+ 77—7rT.i^^l,iV^ + A'7,)^P^Z„ (1) 

Z COS t7w 



'^"'^We use the Pauli metric as discussed by deWit and Smitl 



where 0^ = ta.n~^ (g'/g) is the weak angle, e = gsinOyy is the positron electric charge 
(and hence is positive) and = cos 6'w + sin is the massless photon field. 
Further, = (W^ ^iW^)/2 and Z = —i? sin 6'w + cos 6'w are the massive charged 
and neutral weak boson fields, respectively, while and T~ are the weak isospin 
raising and lowering operators. The vector and axial couplings are 



w 



A' ^ hL{t), (2) 

where tsii'i) is the weak isospin of fermion i ( +i for and z/j; — ^ for di and £j) and 
Qi is the charge of ipi in units of e. 

The three terms in the interaction Lagrangian of Eq.|l| represent the charged- 
current weak interaction, the electromagnetic interaction and the neutral-current weak 
interaction respectively. Note that the Lorentz structure involves only vectors and 
axial vectors; there is no compelling experimental evidence for scalars, pseudoscalars or 
tensor constructions. The minimal Standard Model as described by Eq. |l] successfully! 
explains W and Z decays, neutrino-hadron scattering, neutrino-electron scattering and 
parity-violating electron-hadron neutral-current experiments, providing that radiative 
corrections to order a, the fine-structure constant, are applied. In this chapter we will 
focus our discussion on the charged- current weak interaction in semi-leptonic decays 
involving quarks and leptons from the first family. As an example, consider the decay 
d ue~Ve, for which the T-matrix iJi 

8 + 



where k is the momentum transferred between the hadrons and leptons. For momenta 
small compared to the VT-boson mass, mw, which is typically the case for leptons and 
quarks of the first family, T/j reduces to 



Q 

Tfi = -^nd ^«7m(1 + 75)^d ^e-lA^ + 75)^Fe 



= -^yudJ^Jl (4) 

which is the familiar effective four-fermion interaction with the Fermi constant given 
by Gy/\/2 = (7^/(8m^). Here and are known as the hadron and lepton currents, 
where 



j; = zv.7/.(i + 75)^. = v;: + a^ (s) 



and analogously for J^. 



Table 1. Parameters of the CKM matrix 



Value 




Source 


|Kd| = 


0.9741 ±0.0017 


Nuclear beta decay relative to muon decay 


\Vus\ = 


0.2205 ±0.0018 


f^es and hyperon decays 


\v^\ = 


0.204 ±0.017 


Neutrino production of charm from valence d quarks 


\V,s\ = 


1.01 ±0.18 


D Ke^Ve and D lifetimes 


1 Kfe/Kf, 1 = 0.08 ± 0.02 


Energy spectrum for semileptonic decay of B mesons 


\v,,\ = 


0.040 ± 0.005 


B D*iT7i and B lifetimes 



1.2. Cabibbo-Kobayashi-Maskawa Matrix 

The up-quark, u, obtained in Eq.|^ by the explicit operation of the weak interac- 
tion, 7^(1 ±75), upon the down-quark, d, may not be the physically observed up-quark. 
That is, the quark mass eigenstates are not the same as the weak interaction eigen- 



states The matrix relating the two bases was defined for three families and given an 
explicit parameterization by Kobayashi and MaskawJ. It generalizes the two-family 
case where the matrix is parameterized by a single angle: the Cabibbo ang lei. 

By convention, the three quarks with charge 2/3 (m, c and t) are unmixed and all 
the mixing is expressed in terms of a 3 x 3 unitary matrix V operating on the quarks 
with charge —1/3 {d,s,b): 

( d'\ I Vud Vus Kb \ / \ 

= Vcs V,A \ s\. (6) 

\y j \ Vu Vts Vtb J \ b J 

The values of individual matrix elements can in principle be determined from weak 
decays of the relevant quarks. The Particle Data Group's recommended values! are 
listed in Table |l] except for Vud, which we will have more to say about later in this 
chapter. By construction, the CKM matrix is unitary so its nine elements are not 
independent of each other. Typical parameterizations of the matrix involve three angles 
and one phase. In the Standard Model, the explanation of the observed CP violation 
requires such a phase and hence requires at least three families of quarks. 

The unitarity of the CKM matrix is most stringently tested experimentally in the 
first row of the matrix: 

|Kd|'±|K.|'±|Kfc|' = l. (7) 



"^^In the limit that quark masses are zero, the two bases will be identical and no CKM matrix is 
required. In the lepton sector we will not consider any ramifications of the possibility that neutrino 
masses are non-zero and hence not introduce any CKM matrix in the lepton sector. 



This test is dominated by the accuracy with which Vud can be determined. A failure of 
the test could signal an extension to the Standard Model, such as extra U{1) groups or 
right-handed, 7^(1 — 75), weak- interaction couphngs. We return to these points later. 



1.3. Lagrangian for Nucleons 

In nuclear physics, a similar Lagrangian for nucleons is required with the nu- 
cleons treated as elementary spin-1/2 fermions. The Standard Model interaction La- 
grangian of Eq.[l| may not be immediately applicable for this case, however, since this 
elementary Lagrangian is explicitly constructed for, and tested in, interactions of point- 
like fermions. Nucleons, by contrast, are composite particles of three quarks. To treat 
them as elementary spin-1/2 fermions requires some modification to the elementary 
Lagrangian as can be seen from the following argument. 

Consider the time component of the vector part of the hadron current (see Eq.|^), 
— i^ulii^d = i^ii^d, and evaluate its expectation value between proton and neutron 
states. In the naive quark model, the nucleon comprises three first-generation quarks 
arranged in a fully symmetric way with orbital angular momentum zero and spin 1/2. 
The expectation value is easily evaluatedB: 

{p I V^o' I n) = 1. (8) 

The result of unity implies the vector current normalization in the nucleon is unchanged 
from that in the quark Lagrangian. This, as we shall see, is as it should be and is a con- 
sequence of the conservation of the vector current. In contrast to this, consider the space 
component of the axial- vector part of the hadron current: A'^ = iip^'y-f^ipd — —f^'ipt'^d 
in a non-relativistic limit. Here cr is the Pauli spin matrix, which in its application here 
flips the quark spin. ThuH 

{P I ^3 I n) = -I (9) 

and the normalization of the axial current is indeed modified in a nucleon. The naive 
quark model value of -5/3 is larger than the experimental value of -1.26 obtained from 
the beta decay of the neutron. However, the nonrelativisitic treatment given here is not 
justifiable for quarks as light as 10 MeV. Typical relativistic quark calculations give 
values in the range -1.1 to -1.3 depending on the form of the quark wave functions used 
and the assumed masses of the quarks. 

Since the quark model lacks the predictive precision required, we will follow the 
historic and more phenomenological approach of writing down a general Lagrangian 
using axioms and symmetry principles to limit the number of terms. For example, for 
the weak vector current, V^, the most general form is constructed by our writing down 
all possible ways of making a Lorentz vector out of the vectors at our disposal: p^, the 
incoming momentum of the nucleon; p'^, the outgoing momentum of the nucleon; and 
7^, the spin matrices. If it is further assumed that a nucleon in a nucleus undergoing a 



weak interaction can be treated as a free nucleon, which for the purposes of construct- 
ing interaction operators satisfies the Dirac equation, then the number of terms in the 
current is hmited to three. This latter assumption is known as the impulse approxima- 
tion. Similar arguments also limit the number of terms in the axial current to three. 
Both currents are written as 



and ipn are field operators representing the proton and neutron respectively. The form 
factors g^, gM, gs, S'a, S't, gp are arbitrary functions of a Lorentz scalar constructed 
from the momenta p' and p. It is usual to choose this scalar to be k"^; any other 
scalar can be expressed in terms of fc^ and the nucleon mass, M. The values of these 
form factors in the zero- momentum transfer limit, k"^ —>■ 0, are known as the vector, 
weak magnetism, induced scalar, axial vector, induced tensor and induced pseudoscalar 
coupling constants respectively. 

The number of arbitrary functions in and can be reduced further by requir- 
ing that the currents be Hermitian and invariant under time reversal. These condition^ 
lead to the requirement that gyi^k"^), g^ik"^), g^ik"^) and gp{k'^) be real functions and 
that gsik"^) = 5'T(fc^) = 0. Further, the statement above that the normalization of the 
vector current is unchanged in going from the quark to the nucleon Lagrangian requires 
gvik"^ ^ 0) = 1. Additional constraints on the form factors follow from the conserved 
vector current (CVC) and the partially conserved axial-vector current (PCAC) hy- 
potheses, which we discuss in later sections. 

It is convenient in nuclear physics to make no distinction between protons and 
neutrons, but instead to use an isospin formalism. Let ip = [ipp, ipn) be an eight- 
component spinor operator representing the nucleon field, with ipp and ipn each being 
four-component field operators for the proton and neutron respectively. Further, intro- 
duce Pauli spin matrices, r^, and r^, forming the Cartesian components of a vector, 
T, that are generalized in block form to act on ip. For example, acting with on a 
pure proton state or a pure neutron state leads to 





(11) 



These expressions serve to specify our isospin conventions. The currents in Eq. |10] in 
an isospin formulation are now written 



v: 



h 




2M 



4 = z^[^^(P)7^75 + ^|^a^.A;,75 + ^^?p(A:')A:^75]^r>, (12) 

where stands for ± ir^ with a plus sign for n — > pe~Vf. and a minus sign for 
p ne^Ve charge-changing transitions. 

Another symmetry of these currents to be investigated is their behaviour under 
charge conjugation, which interchanges particles and antiparticles, coupled with a 180° 
rotation about the 2-axis in isospin soace. This combination is called G-parity and is 
a symmetry of the strong interactions^: 

G = Cexp(«7rT2). (13) 
Now, under charge conjugation, we havei 

cv^.ct = ^ ^,4, (14) 

where i = p,n and the representation of C in terms of the Dirac matrix, 72, is specific 
to the metric we are using. Further, under an isospin rotation. 

Thus, combining these two operations, we have 

GijpG^ = -ij^ ; GiJnG^ = i:^. (16) 
Applying these steps to the first term in the vector current, Eq. [TO, we get 



= -«(01(747m)"/3(V^p)/3 
= +^(^p)M747M)?a(€)^ 

= +#j727j7r72^n 

= (17) 

where, in the third line, the expression is written out in component form and reordered 
in the fourth line, introducing a minus sign, because the fermion field operators anti- 
commute. Thus, the first term in the vector current is invariant under G-parity trans- 
formations. Applying the same steps to the first term of the axial-current in Eq. ID, we 
get 



GAlG^ = lG^^G^l^l,G^r^G^ 

= -4- (18) 

We considered iust the first terms of Vj!' and A^, because tliese are tlie only terms in 
tlie quark currents of the Standard ModeL WeinbergtU has termed currents with the 
standard quark-model transformation properties under G-parity as "first class" , while 
currents with the opposite properties, 



QyHn)Q\ ^ _yh(n) 



^^j(//)G't = Al^^^\ (19) 

are termed "second class". Under this classification, the vector and weak magnetism 
terms of Vj^ and the axial-vector and pseudoscalar terms of A^^ are first-class currents, 
while the induced scalar term of V!^ and the induced tensor term of A^, are second- 
class currents. To date, there is no unambiguous evidence for the presence of second- 
class currents in nuclear beta decay, but it is nonetheless an interesting experimental 
challenge to place limits on the possible presence of these second-class currents. 

2. Vector Interaction 

2.1. Conserved Vector Current Hypothesis 

Briefly, let us return to the quark Lagrangian of Eq. |l] for the electromagnetic 
interaction and, in analogy with Eq. ^, write down the electromagnetic current for the 
first family of quarks, 

vr = ^¥ui,i^u - (20) 

in units of the coupling constant, e, which is the positron electric charge. Again, we 
can introduce isospin by defining a quark doublet spinor ipq = {ipu,ipd) such that 

V^"' = iQ+ii,i,i^,I + *g_^,7^V^,r3, (21) 

where / is a unit matrix in isospin space; Q+ is the average quark charge, Q+ = ^(| — ^), 
and Q- the average difference, Q- = ^(| -|- |). Note that the electromagnetic current 



has both isoscalar and isovector components. When the steps outlined in sect. |1.3| are 
repeated, the electromagnetic current between nucleons likewise exhibits isoscalar and 
isovector components and is written 



irem 



2M 
Flip) 
2M 



(22) 



where now is the average nucleon charge and Q- the average charge difference: 
(5+ = Q- = ^. Again, symmetry principles hmit the number of terms to three for both 
the isoscalar and isovector currents. In terms of proton and neutron form factors, we 
identify 



FfiP) = Ff{e) + FJ^{k% 

Fne) = F[{k') - Fpik'), z = 1,2,3. (23) 

Further, in going from quark to nucleon Lagrangians, the normalization of the vec- 
tor current remains unchanged so that the first form factor Fi in the fc^ — limit 
simply specifies the charge of the nucleon in units of e: viz. -Ff (0) = 1, F^{0) = 0, 
or equivalently F^{0) = 1, F^{0) = 1. The second form factor relates to the anoma- 
lous magnetic moment of the nucleon (in units of the nuclear magneton, e/(2M), with 
h = c = 1): viz. -Ff (0) = fip — 1 = 1.79 and F2{0) = fin = —1.91, or equivalently 
^2^(0) = —0.12 and -^2^(0) = 3.70. Finally, the electromagnetic current, V^^*" is a con- 
served current: that is, c^^V^™ = or, as written in momentum space, fc^V^*" = 0. 
This follows from the fact that the electromagnetic Lagrangian is invariant under the 
phase transformation ijj ^ ip' = e^^^tp and ip ^ ip = e'^'^^ip, where A is a real con- 
stant. From Noether's theoremi, for every transformation that leaves the Lagrangian 
invariant there is a corresponding conserved current, in this case the electromagnetic 
current. Evaluating k^V^"^ term by term we have: 

W)K^Mp) = WhM - P'M^P) = iW)iM - M)i,{p) = (24) 

from the Dirac equation of a free nucleon; and cr^uk^k^ summed over repeated indices 
is zero from the antisymmetry of cr^^; while k"^ = k"^ ^ 0. Thus, current conservation 
requires that 

F3{k^) = (25) 

for both isoscalar and isovector currents. 

It is noted that the isovector electromagnetic current, Eq. ^ and the weak vector 
current, Eq. 0, are identical in structure with the isospin matrix in the former case 
replaced by the matrix r" in the latter case. Feynman and Gell-Mann0 postulated that 
these currents are members of an isotriplet vector of current operators. A consequence 
of this hypothesis is that the weak vector current is also a conserved current. This is 
called the conserved vector current (CVC) hypothesis, from which it follows that 



g^(e) = Fl^ik^) 1 as e ^0 
g,,{e) = F^ik') 

gsiP) = F^ik') = 0. (26) 



Note that gv{0) = 1 always, irrespective of the nucleus in which it might be experi- 
mentally measured. That is, the vector coupling constant is not renormalized in the 
nuclear medium. Likewise, the weak- magnet ism form factor is uniquely related to the 
corresponding electromagnetic form factor: for beta decays between isospin analogues 
it relates to the anomalous isovector magnetic moment of the nuclear states in ques- 
tion, while for non-analogue transitions it relates to the corresponding isovector Ml 
gamma-decay rate. 

So far, experimental tests of CVC have taken two forms: (^zj showing 5'v(0) is a con- 
stant irrespective of the nucleus in which it is measured, and (ii) showing gM^O) = -^2^(0) 
by measuring gM{0) from either the lepton energy spectrum or an angular correlation in 
a beta-decay experiment and comparing it with -^2^(0) obtained from the lifetime and 
branching ratio for the analogous 7-transition. We discuss both tests in the following 
sections. 

2.2. Superallowed Fermi Transitions 

Precise measurements of the /5-decay between nuclear states with [J'",T) = 
(0^, 1) provide the most straightforward yet demanding probe of the vector interaction 
in nuclei. Because the axial current cannot contribute in first order to transitions be- 
tween spin-0 states, the experimental /t-value is related directly to the vector coupling 
constant, = GpVud- 

ft(l + 5n) = o, (27) 

with 



(Mv)^ = 2(1 - 6c) for T = 1 
K/{nc)^ = 2'n'Hn2h/{meC^f = (8120.271 ± 0.012) x 10"^°GeV-^s, (28) 

where the physical constants used to evaluate K were taken from ref. i. Here / is the 
statistical rate function, t is the partial half-life for the transition, and (My) is the Fermi 
matrix element. The radiative corrections, 6^ and A^, and the charge correction, Sc, are 
each of order 0.01, with and particularly 6c being dependent on nuclear structure. 
There are theoretical uncertainties introduced via these correction terms but they are at 
least another order of magnitude smaller, so Eq. ^ provides an experimental technique 
for determining Gy to better than a part in a thousand. 

To date, nine 0"^ —>■ 0+ transitions have been studied to this level of precision; 
they cover a wide range of nuclear masses from ^°C, the lightest parent, to ^^Co, the 
heaviest. The results have been used to test the constancy of G^, which is a prediction 
of CVC, and to derive the value of Vud, the up-down element in the Cabibbo-Kobayashi- 
Maskawa quark-mixing matrix. A unitarity test on the matrix, made possible by these 
results, is a fundamental measure of the adequacy of the three-generation Standard 
Model. 



Table 2. Experimental results {QeCi ^i/2 branching ratio, R) and calculated corrections 
{Pec-, and 5r) for 0+ — > 0+ transitions. 





WEC 


tl/2 


it 


^EC 




oc 


Or. 






(keV) 


(ms) 


(%) 


(%) 


(s) 


(%) 


(%) 


(s) 


IOC 


1907.77(9) 


19290(12) 


1.4638(22) 


0.296 


3040.1(51) 


0.18(4) 


1.30(4) 


3074.0(54) 


14Q 


2830.51(22) 


70603(18) 


99.336(10) 


0.087 


3038.1(18) 


0.23(4) 


1.26(5) 


3069.2(28) 


26m 


4232.42(35) 


6344.9(19) 


> 99.97 


0.083 


3035.8(17) 


0.31(4) 


1.45(2) 


3070.2(23) 


34C1 


5491.71(22) 


1525.76(88) 


> 99.988 


0.078 


3048.4(19) 


0.54(4) 


1.33(3) 


3072.1(25) 


38m 


6043.76(56) 


923.95(64) 


> 99.998 


0.082 


3047.9(26) 


0.56(4) 


1.33(4) 


3071.1(32) 


42SC 


6425.58(28) 


680.72(26) 


99.9941(14) 


0.095 


3045.1(14) 


0.36(4) 


1.47(5) 


3078.7(24) 


46y 


7050.63(69) 


422.51(11) 


99.9848(13) 


0.096 


3044.6(18) 


0.37(4) 


1.40(6) 


3075.8(28) 


50Mn 


7632.39(28) 


283.25(14) 


99.942(3) 


0.100 


3043.7(16) 


0.41(4) 


1.40(7) 


3073.8(28) 


^4 Co 


8242.56(28) 


193.270(63) 


99.9955(6) 


0.104 


3045.8(11) 


0.50(4) 


1.39(7) 


3073.1(28) 



Average, J't 3073.1(11) 
xV(iV-l) 1.34 



2.2.1. Experiments 

The ft value that characterizes a particular /5-transition is determined by three mea- 
sured parameters: the transition energy, Qec: which is used in calculating the statistical 
rate function, /; the half-life, ti/2, of the /9-emitter and the branching ratio, R, for the 
transition of interest, which together yield the partial half-life, t. World data have been 
surveyed recentlyB for eight of the best-known superallowed emitters: ^^O, ^^'"Al, ^^Cl, 
38mj^^ ^^Sc, ^^V, ^°Mn and ^^Co. The resul ts giv en in Table 2 are drawn from that sur- 
vey, augmented by data on decay — — iMi'S'El'Eli^ill and by several other recent 



measurements of significance .B'B'EJ'EllHSEZ 

The experimental uncertainties quoted for the measured decay energies (typically 
1 part in 20,000), half-lives (1 part in 3,000) and branching ratios (1 part in 70,000) are 
at the limits of what is currently possible. Even though such quantities are frequently 
measured in nuclear physics, the extreme demands on precision in these studies have 
led to experimental techniques highly refined for the purpose. Several examples, briefly 
described, will illustrate the point. 

2.2.1.1 Decay energies. The nine superallowed /3-emitters listed in Table 2 popu- 
late daughters with stable ground states. As a result, reaction Q- value measurements, 
which lend themselves to experimental precision, can be used to determine decay ener- 
gies instead of direct measurements of positron end-point energies, which are frought 
with precision-limiting difficulties. For example, in the case of 26^^! — i>^^Mg, the decay 
energy, Qeci has been measured by three different routes: 

(i) ^^Mg(p, n)^^'"Al. The threshold energy for this reaction was determinedil by 
detection of the ^sm^^j activity produced as a function of proton energy, with a 26 Mg 
target located in the focal plane of a magnetic spectrometer. The energy of the proton 



beam was calibrated by comparison with a ^'^'^Cs"*" beam accelerated through a voltage 
difference V to pass along the same path in the spectrometer. The value of V was then 
compared via a resistive divider with a one- volt standard. The measured Q value, Qpn, 
derived from the observed threshold energy, is related to Qec by: 

Qec = -Qpn - 782.354 keV. (29) 

(ii) ^^Mg(p, 7)^^™'A1 and ^^Mg(n, 7)^^Mg. For the first reaction, proton ener- 
gies at four resonances were establishedc3 by comparison with known resonances in 
^^Al(p, 7)^^Si. For the second, thermal neutrons were used, with 7-ray energies being 
determinedc3 by comparison with known calibration lines. The measured Q values yield 
Qec through the relation: 



Qec = Qn^-Qiry- 782.354 keV. (30) 

Though comparable uncertainties are claimed as for method (i), it is evident that this 
method is less direct and involves a variety of secondary calibration standards. 

(Hi) ^^Mg('^ife, t)^^Al (1058 keV). This reaction was observed^ in conjunction 
with ^^N(^ife, t)^^0 from a composite target in which ^''N had been implanted into 
^^Mg. The target assembly was biased periodically to ~80 kV, thus permitting the two 
triton groups corresponding to the ^^O ground state and the 1058 keV state of ^^Al 
alternately to follow the same path through a magnetic spectrometer. The voltage at 
which they perfectly overlapped corresponds to the Q-value difference, AQ, between 
the two reactions. The energy difference between the 1058 keV state and the ^^Al 0"*" 
isomer was determined^ from the relevant 7-ray energy, E^-^. Similarly, the ground 
state of ^''N can be related to the 0~^(T = 1) excited state by a 7 ray of energy, E^^. 
Thus, 



Qsc(A = 26) - Qec{^ = 14) = E^, - E^, - AQ. (31) 

This method, too, depends on secondary calibration standards but it provides valuable 
interconnections betwen pairs of superallowed 0"*" —>■ 0^ transition Qec values. 

2.2.1.2 Half-lives. To achieve the required precision, isotopic purity of the sample 
must be assured. Currently, all emitters except ^^O and ^''C have been measured^ 
with an on-line isotope separator. Typically, thousands of multi-scaled spectra are 
acquired from successive samples observed in a large-solid-angle, high-efficiency, low- 
background gas proportional counter; the counting electronics have a well-defined, non- 
extendable dead time. The best analysis procedure uses modified Poisson variances and 
is based on a simultaneous fit to up to 500 decay curves with the same half-life but with 
individual intensity parameters for each decay curve. This procedure is checked^ with 
hypothetical data, computer-generated to simulate closely the experimental counting 
conditions, and has been proven correct at the 0.01% level. 

2.2.1.3 Branching ratios. In most cases, the superallowed transition dominates 
(> 99% ) the decay of the parent nucleus and is best determined via measurement 



of the other decay branches, their branching ratios being subtracted from 100% . For 
these other branches, precision is not as important as completeness; therefore special 
techniques have been developed {e.g. ref.i§) to observe /3-delayed 7 rays from branches 
as weak as a few parts per million. For the only case in which the superallowed transition 
is itself weak - the decay of ^''C - the technique usedlll is to compare the intensity of 
/3-delayed 7 rays in ^°B, produced via the reaction ^''B(p, n)^°C (/9"'")^°B*, with the 
same 7 rays produced via ^°B(p, p')^°B*; for added precision, a large 7-ray arrayill has 
now been used to observe these 7 rays in alternating off-line, on-line measurements. 

2.2.2. The ft values 

The statistical rate function, /, is defined by 

f = pW{Wo-WfF{Z,W)C{W)dW, (32) 

where W and p are the electron energy and momentum, respectively, in electron rest- 
mass units, and Wq is the maximum /^-energy, which is related directly to the measured 
Qec via the equation Wo + 1 = QEc/{.i^e(?)- The quantity W) is known as the 
Fermi function and C{W) the shape correction factoiil. The latter has a value close 
to unity for most values of W , its departure from unity reflecting the influence of 
screening by atomic electrons, the dependence of the nuclear matrix elements on W, 
and the presence of second-forbidden matrix elements. Values of the statistical rate 
function have been computed for each superallowed transition following the techniques 
described in refs.Bil. Each result has an associated uncertainty reflecting the experi- 
mental uncertainty in the Qec value used in the calculation. 

Next, the partial half-life for each superallowed emitter must be obtained from 
its total half-life and branching ratio according to the formula 

t = h,2{l + PEc)/R. (33) 

where Pec is the calculated electron-capture probability. Table 2 lists the Pec values 
calculated by the method described in ref.lli, which employs exact bound-state elec- 
tron wave functions. Though the calculated values of Pec are small in all cases, the 
correction is comparable to the experimental uncertainties on ti/2 and R, and must be 
properly accounted for. 

Finally, the /t-values that incorporate all experimental uncertainties are listed 
for the nine emitters in Table 2. They incorporate the complete set of world data and 
form the experimental basis for the fundamental tests that follow. 

2.2.3. Theoretical corrections 

In a nucleus, the charged particles involved in (3 decay interact with an external electro- 
magnetic field, which leads to effects of order (Za)"*, m=l,2..., where Z is the atomic 



number of the daughter nucleus and a is the fine-structure constant. For the electron 
(or positron), these effects are accounted for in the calculation of / (see Eq.^), which 
usesliii exact solutions of the Dirac equation for electrons in the field of a modelled 
nuclear charge distribution. For the proton, these effects are reflected in the breaking of 
nuclear isospin symmetry and a resulting correction, be-, to the Fermi matrix element. 

There are also effects that arise from interactions between the charged particles 
themselves. These are of order Z^a", with m < n, and lead to the radiative corrections, 
Sr and A^. By convention, the former incorporates those corrections that are nuclear- 
structure dependent, the latter, those that are not. 

2.2.3.1 Charge correction, 5c- Both Coulomb and charge- dependent nuclear forces 
destroy isospin symmetry between the initial and flnal states in superallowed /3-decay. 
The consequences are twofold: there are different degrees of conflguration mixing in 
the two states, and, because their binding energies are not identical, their radial wave 
functions differ. CalculationffliS accommodate both effects by splitting 5c into 
two components, and 5c2 , to account for configuration mixing and radial mismatch, 
respectively. In all cases 5c2 > • 

To establish 5ci, shell- model wave functions were calculated with appropriate 
nucleon-nucleon interactions modified to include: (i) two-body Coulomb terms in the 
proton-proton part of the Hamiltonian, (ii) a 2% increase in proton-neutron elements 
(justified by the charge dependence observed in nucleon-nucleon scattering data), and 
(Hi) one-body elements determined from the single-particle energies of closed-core- 
plus-proton and -neutron nuclei. These basic elements were then adjusted to reproduce 
the b- and c-coefficients of the isobaric-multiplet mass equation that describes the 
multiplet involved in the superallowed transition. The results have also been tested 
independently by measurements^, in some of the same nuclei, of weak non-analogue 
0+ O"*" transitions, which are sensitive to the same effects. 

The larger of the two components, 5c2, has been calculated from two different 
approaches. OneH describes the proton in the parent nucleus and the neutron in the 
daughter nucleus by full parentage expansions in terms of Woods-Saxon wave func- 
tions suitably matched to the experimentally known separation energies. The other0 
employs a self-consistent Hartree-Fock calculation in which the overall strength of the 
central potential was adjusted so that the single-particle eigenvalues agree with the 
experimental separation energies. The former approach leads to values for 5c2 that are 
consistently larger than those from the latter. 

In Table 2^we list 5c {=5ci +^€2) for ^^^^ transitions. The values for 5c\ were 
taken from ref.lll updated by new results in ref.H; those for 5c2 are the unweighted 
averages from the results of refs. IS. "Statistical" uncertainties, reflecting the scatter 
of the calculated values, were derived in ref. and have been adopted here. 

2.2.3.2 Radiative corrections, 5r and A^. In the context of Standard-Model tests, 
we need consider only radiative-correction terms that distinguish semi-leptonic nuclear 
/3-decay from purely leptonic /i-decay0. By grouping them according to whether they 
are nuclear-structure dependent {5r) or not (A^), we writeS: 



a 

5r = —[9iEm)+62 + 6s + 2CNs] (34) 

ZTT 

a 

Al = — [41n(mz/mp)+ln(mp/mA) + 2CBorn] + ..., (35) 

ZTT 

where, for clarity, only leading-order terms have been presented for A^. In these equa- 
tions, Em is the maximum electron energy in the /5-decay, rrip and rriz are the masses 
of the proton and Z-boson, and is a mass parameter associated with the Ai-meson 
but, in this application, is allowed the range 400 < < 1600 MeV. The func- 
tion g{E, Efn) was first defined by Sirlinil as part of the order-a universal photonic 
contribution arising from the weak vector current; it is here averaged over the elec- 
tron spectrum to give^(-Em). The effects of order-Za^ and -Z'^a^ have recently been 
definitively calculated^S and are incorporated as 62 and ^3, respectively. Finally, the 
terms CBorn and C^s arise from the order-a axial- vector photonic contributions: the 
former accounts for single-nucleon contributions, while the latter covers two-nucleon 
contributions and is consequently dependent on nuclear structure. We adopt here the 
calculations of Towner^. 

Built from these components (see Eq. ^4]), calculated values for Sr are listed in 
Table 2. The largest component, 'g{Em), is responsible for approximately 1% in each 
case, S2 and ^3 together account for less than 0.5%, and the effects of Cns range from 
-0.3% to +0.1%. The calculated value for A^ is 



A^ = (2.46 ± 0.09) %. (36) 

In this case, the leading logarithic term dominates, with a contribution of 2.1%. The 
remainder arises from CBorn (0.20%) and higher order terms not explicitly listed in 
Eq. ^ The second term in the equation effectively establishes the quoted uncertainty, 
which reflects the range of values allowed for m^. 



2.2.4. Implications for the Standard Model 

It is now convenient to define a "corrected" J-'t value, which is related to the conven- 
tional ft value via: 

J't = ft{l + 6R){l-6c). (37) 

As a comparison with Eg . |27| makes clear, all nuclear-structure-dependent effects should 
have been removed from the J-'t values. It is the consistency of the determined J-'t values 
over a wide variety of superallowed O"*" — ^ 0^ transitions that constitutes a demanding 
test of CVC. 

Table 2 contains all requisite experimental and calculated information for the 
determination of J^t values, and the results for all nine transitions appear in the last 
column of that table together with their weighted average, J-'t, and scaled "statistical" 
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Fig. 1. J^t values for the nine precision data and the best least squares one-parameter fit 



uncertainty. The same results are displayed in Fig.|l|. There is no statistically significant 
evidence of inconsistencies in the data, thus verifying the expectation of CVC at the 
level of 4 X 10~^, the fractional uncertainty quoted on the average jFt-value. 

With the consistency of the measured jFt- values established, we turn to the eval- 
uation of Vud and the unitarity test of the CKM matrix. The up-down quark mixing 
matrix element is given by 




_ K 

2Glil + Al)W ^^^^ 

In using the average J-'t value from Table 2 in Eq.^ it is important now for us to 
include in the jFt-value uncertainty the effects of a "systematic" uncertainty in the 
charge correction, 6c, as well as the "statistical" uncertainty already described. The 
"systematic" uncertainty of ±0.08% arises from systematic differences between two 
independent model calculationJUa of 6c2, and has been elaborated on in ref. El. The 
result is 

JY = 3073.1 ±3.0 s (39) 

Thus, with 



G^/iUcf = (1.16639 ± 0.00002) x 10"^ GeV"^ (40) 
from ref. i and from Eq. we obtain 

Vud = 0.9736 ± 0.0006. (41) 

The quoted uncertainty is dominated by uncertainties in the theoretical corrections, 
A^ and 6c', experimental errors on the input data for the jFt- value contribute less than 



±0.0002 to the uncertainty in Eq.[41 



This result for Vud leads to the most demanding test yet available of the unitarity 
of the Cabibbo-Kobayashi-Maskawa matrix, in which 



\V^d\^+\Vus\^+\Vub\'=l (42) 

must be satisfied if there are three generations of quarks prescribed by the Standard 
Model. Evidently, Vud is by far the largest of the three matrix elements and its uncer- 
tainty the most critical. The value of Vus is taken to bell Vus \ = 0.2205 ± 0.0018, the 
average of two results, one from the analysis of K^s decays, the other from hyperon 
decays. Finally, V^b can be derived, with substantial uncertainty, from the semi-leptonic 
decay of B mesons, but is sufficiently smalli, at 0.0032 ± 0.0009, as to be negligible in 
the unitarity test. The result. 



I Vud I '± I Vus \^+\Vub\' = 0.9965 ± 0.0015, (43) 

differs from unitarity at the 98% confidence level. 

The significance of this apparent non-unitarity is still being debated. Taken at 
face value it indicates the need for some extension to the three-generation Standard 
Model (see sect.^). However, it may also refiect some undiagnosed inadequacy in the 
evaluation of | Vus | ^ or | Vud I • For example, it has been proposed^! that the nuclear- 
model-dependent calculations of 6c should be replaced by an empirical approach in 
which the uncorrected data are fitted by a smooth Z-dependent function. Extrapolation 
of this function to Z = would yield the recommended value of J-'t. With the data 
in Table 2, a variation of this approach is to allow a Z-dependent fit of the nine J-'t- 
values: doing so with linear Z dependence does lower the reduced somewhat (from 
1.34 to 1.04). The result is J-'t{Z = 0) = 3068.1 ±3.1, which yields a new unitarity 
sum (Eq. ^) of 0.9980 ± 0.0015. It must be emphasized, however, that although this 
approach does yield a less provocative result, there is no substantive explanation for 
any correction not already accounted for in 6c and 6r, and there is no statistically 
significant indication from the data that one is required. 

'^'^A larger value is recommended in the analysis of ref@ 



2.3. Experimental Tests of CVC in Weak Magnetism Terms 



The conserved vector current hypothesis, besides requiring that QyIO) is constant 
(= 1), as verified in superallowed transitions, also predicts that the weak magnetism 
form factor is uniquely related to the corresponding electromagnetic form factor. In 
this section we will discuss experimental tests of this latter prediction, concentrating 
our attention principally on the {J'^, T = l'^, 1) (0"'', 0) allowed axial-vector mirror 
/5-transitions in the A = 12 nuclei: ^^B ^^C +e~ +Ve and ^^N ^^C +e+ + u^. For 
oriented parent nuclei, the beta-decay rates (upper sign for ^^N decay and lower sign 
for ^^B decay) are given byS 

d'T{e^) = l.^v;^,\F^{0)\' F±iZ,W)pW{Wo-W)\l + v± + a±W)dWdn, 

x[l ± (/ii - /i_i)(l + a±iy)Pi(cos^e) + (1 -3/io)a±iyP2(cos^e)],(44) 

where F±{Z,W) is the Fermi function for the decays, p and W are the electron 
momentum and total energy in electron rest-mass units; and Wq is the maximum 
electron energy. The shape correction factor, C{W), has been explicitly written as 
(1 + '7± + a±W). In addition, hi, h-i, ho are the populations of the m = 1, —1, states 
of ^^N or ^^B, normalized such that hi + + /iq = 1; cos6'e = p ■ z, and F = ln2/t, 
where t is the decay half-life. The coefficients 7]±, a± and a± are expressed, in the 
formalism of Hwang and Primakoffii, in terms of nuclear form factors, F: 



V± = 



W, 



■± 



a± = =F 




SMFfiO)^^^ ^ 



1 (Fi,{0) F|(0)\ 

where M is the nucleon mass and 6±{aZ) are radiative corrections calculated to be 
5_ = 1.31 GeV"^ and 6+ = 2.69 GeV"^ by Behrens and SzybiszB. These nuclear form 
factors originate in an elementary-particle treatment^ of the decays in which no explicit 
consideration is given to the dependence of the wavefunctions on the coordinates of the 
constituent nucleons, nor is the weak current itself decomposed into contributions from 
individual nucleons. Rather, the details of nuclear structure are buried in the nuclear 
form factors, which have then to be determined from experiment. 

With sufficient data, this approach provides a treatment that is essentially inde- 
pendent of nuclear models and permits analyses to explore the consequences of sym- 
metry properties and conservation hypotheses. But with insufficient data, some of the 
form factors have to be obtained from other sources, the most common being to appeal 



to the impulse approximation. In the latter case, the starting point is the single-nucleon 



currents of Eq. |T2| reduced to their nonrelativistic form to yield one-body operators, the 
matrix elements of which are evaluated with many-body wavefunctions. The coupling 
constants of Eq.|T2|, g{k'^) ~ g{0), are determined from data on the free proton and 
neutron and treated as constants, unchanged from nucleus to nucleus. The connection 
with the nuclear form factors of the elementary-particle approach is given bySJ^l 

v/2F±(0) ^ -(1 + (7m)('t)-(L) 

y2F|(0) ^ -{g^±g^){a)~2gj,{zr{(T-p)), (46) 

where, for example, (cr) is the Gamow- Teller matrix element defined more completely 
below in Eq. 0. 

The transition rate for the analogous gamma decay, ^^C*(l^,l) ^^^C(0^,0), is 
specified by 

r^c*- ''C) = |^|V2M0)P, (47) 

where a is the fine-structure constant, the photon decay energy and /i(0) the 
isovector magnetic nuclear form facipr. which, in the nonrelativistic limit, connects 
to impulse-approximation quant it ieseS'Ea via 

V2/i(0)^-(/ip-/i„+l)(o-)-(L), (48) 

with Hp and /i„ the proton and neutron nucleon magnetic moments respectively. The 
test of CVC consists of the verification of the relation 



F±(0) = v^^(O). (49) 

The y/2 displayed in Eq.|9| occurs because the isospin operator for beta decay is (r^ ± 
ir^), whereas the isospin spherical tensor operator that relates to the T3 operator of 
gamma decay through an isospin rotation is (r^ ± ir^)/-\/2. 

The weighted average of six measurements^ of that use the technique of 
nuclear resonance fluorescence (after reevaluation with currently accepted branching 
ratios) is = 40.4 ± 2.0 eV. The three most recent measurements of F^ that use 
^^C(^e') yield F^ = 37.0 ± 1.1 eV (ref.@), 35.74 ± 0.86 eV (ref.i) and 38.5 ± 0.8 eV 
(ref.li). Following de BraeckeleerEl, we adopt the value F-y = 38.8 ± 1.6 eV, which is an 
average of the combined data with an error chosen to make the two data sets consistent 
with each other. There is a small correction!^ due to isospin mixing of the ^^C*(l"^, 1) 
state at 15.11 MeV with the (1~'',0) state at 12.71 MeV, but its magnitude is smaller 



than the stated error on so we consider it to be included in that estimate. Thus, 
from Eqs.^ and ^ we obtain as a prediction of CVC: 



F^(0) = 2.02 ±0.04. (50) 

Note that the form factors Fa and Fe both contain second-class current terms 
through the impulse-approximation coupling constant, g^. Thus, the analysis proceeds 
in one of two ways: either one accepts the CVC prediction, Eq. |0], and uses additional 
data to set limits on second-class currents, or one sets = and uses additional data 
to verify the CVC hypothesis. Following^ the latter choice, and deducing axial form 
factors from the experimental half-livesEll to be -^4(0) = 0.512 ± 0.001 and -^4(0) = 
0.481 ± 0.002, the slope parameters, a±, (see Eq.^) are predicted from CVC to be 

a_ = 5.61 + 1.31 = 6.92 GeV"^ 
a+ = -5.97 + 2.69 = -3.28 GeV"^ 

Aa = a_-a+ = 10.20 ± 0.16 GeV"\ (51) 

The order-aZ correctional, 6±{aZ), which appear explicitly as the second terms in the 
equations for a±, are significant and shift the prediction away from the symmetric result 
of a_ ~ — a_|_ expected in the absence of any isospin-symmetry breaking. Behrens and 
SyzbiszB also predict a sizeable quadratic term in the shape correction factor, b±W'^, 
in Eq. which should be included in the analysis. However the contribution from this 
term is roughly the same in e~ and e"*" transitions, so it can be omitted providing we 
restrict our attention to Aa, which is more reliable experimentally and theoretically 
since many systematic uncertainties cancel out. 

The CVC hypothesis alone cannot predict values for the correlation coefficients, 
a±, because of the presence of the Fe form factor in Eq. However a prediction can 
be given for the difference!^: 

Aa = a_ - a+ = 2.89 ± 0.04 GeV"\ (52) 

where, as in Eq.|^, the quoted error reflects only the uncertainty in the Ml 7-decay 
width, r^. The experimental test of the CVC hypothesis, therefore, is to measure the 
difference in slope parameters, Aa, or the correlation coefficients, Aa, in the mirror 
mass-12 transitions and compare with the predictions in Eqs.|51| and We consider 
some examples in the following sections. 

2.3.1. The total decay rate 

For non-aligned ^^B and ^^N nuclei, we can integrate over electron directions and en- 
ergies in Eq.^, giving for the ratio of decay rates in the mirror transitions the result 



/t+ ^ |Fi(0)pl+r7, + a„iyo-/2 
ft- |F+(0)|n + r/+ + a+iyoV2 

l-^^^^±i^^ + (53) 

where / is the integrated Fermi function (see Eq.|3^). The energy- dependent term, aW, 
has been replaced by its average value, aWo/2, and terms proportional to (W(^ — Wq) 
have been neglected. The second line represents the impulse-approximation limit. It 
would seem that a comparison of ft values in mirror transitions would be a profitable 
way to search for second-class currents. This, indeed, was proposed by Wilkinsonil and 
first results appeared encouraging. However, a definitive experiment was performed by 
Wilkinson and Alburgeril in which the asymmetry of ft values was measured as a 
function of {W^ + Wq) in the beta decays of ^Li and ^B. This is a unique case in which 
the daughter state is broad, decaying itself by a emission. This effectively permits the 
beta end-point energy to be varied within a single decay. Their results showed the ratio 
of ft values to be almost independent of (Wq' + Wq), implying little or no evidence 
for second-class currents. 

Another difficulty with this type of analysis concerns the existence of the addi- 
tional term 5nuci in Eq.|5^. It is given by | (cr)~/(cr)+ p= 1 + and arises from the 
fact that (cr)"*" ^ (cr)~ since the wavefunctions for -^^B and ^^N are not exact isospin 
eigenstates: the isospin symmetry is broken by the internucleon electromagnetic interac- 
tion and the neutron-proton mass difference. Nuclear-structure calculationJ^El indeed 
compute values for 5nuci that are comparable with the total experimental asymmetry, 
and thus mask the evidence, if any, for the presence of second-class currents. 



2.3.2. Electron-spectrum shape 

To study the electron-energy spectrum for non-aligned ^^B and ^^N nuclei, we integrate 
only over electron directions in Eq. (|4^) . A measurement of the spectrum makes it 
possible to determine the coefficient of the energy-dependent term a±W. This is the 
classic experiment of Lee et al. which first claimed to verify CVC However, a re- 
appraisal by Calaprice and Holsteinii questioned the original analysis of the experiment 
and weakened the experimental support for the validity of CVC. The critique has been 
discussed in a new analysis by the original group of experimentersil using the data 
from the original experiment, with the conclusion that their data still support CVC. 
At the same time, a new measurement of a± was performed by Kaina et a/. S using 
a Nal crystal as a spectrum analyzer instead of the magnetic spectrometer employed 
in the earlier measurements. The results for a± were not in agreement with previous 
(revised or not) determinations. 

One cause of the trouble lies in the unresolved decay branch from ^^N to the first 
excited state in ^^C. The conclusion derived from the experiment depends strongly upon 



Table 3. Difference of slope parameters, Aa, for B and N mirror beta decays 

Aa bi Remarks 
( GeV-^) (%) 

9.2 ± 1.2 2.4 ref.0, as revised by0 and a+ increased by 0.7 

15.3±2.8 2.4 ref.iil, as revised by0 and a+ increased by 0.7 

8.6 ±2.4 2.1 ref.§ 

9.1 ± 0.9 1.8 ref.B with a+ decreased by 0.7 

10.5 ± 2.1 ref.l^, branching ratio, bf, not given 

"Branching ratio assumed in original reference. 



the knowledge of the branching ratio, bf, for the decay to that state, which must be 
subtracted from the observed result. For example!^, a change from 1.8% to 2.1% in bf 
produces a variation in the extracted a+ of —0.7 GeV^^, a 25% change. A second cause 
of trouble lies in the systematic uncertainties surrounding the experimental energy- 
dependent response functions. So far, these seem insurmountable in the determination 
of a±, but are deemed largely to cancel in Aa; as a result, it is the latter quantity 
that we quote here. In Table |^ are the results from four experiments as corrected for 
a common 2.1% branching ratio as reviewed by GrenacsS, and a preliminary resultii 
from a fifth experiment. The weighted average is Aa = 9.5 ± 0.7. By combining in 
quadrature the effects of an uncertainty of 0.3% in bf, we produce a recommended 
value of 

Aa = 9.5 ± 1.0 GeV-\ (54) 
which agrees with the CVC prediction given in Eq.|5l|. 

2.3.3. Polarization and alignment measurements 

We now turn to experiments with oriented nuclei and measurements of the decay 
electron's angular distribution as given in Eq. 



W{9e) = l±ihi- /i_i)(l + a±W)Pi{cos9e) + (1 - 3ho)a±WP2icos9e). (55) 

Typical experimental arrangements are described by MinamisonoH, Haskell and Mad- 
anskyS and Haskell et a/.B Polarized ^^B and ^^N recoils are produced through the 
reactions ^°B(^He, n)^^N and ^^B((i,p)^^B, with the precise magnetic substate popu- 
lations of the product nuclei depending on the incident beam energy and the recoil 
angle. The polarized recoiling nuclei are implanted into monocrystalline metal foils in 



Table 4. Beta-decay asymmetry coefficients, a± in GeV 



a+(^^N) Reference 



-(0.07 ±0.20) 
+ (0.24 ±0.44) 
+ (0.25 + 0.34) 
+ (0.10 + 0.30) 
-(0.01 + 0.21) 
+ (0.046 + 0.053) 
-(0.174 + 0.059) 

-(0.043 + 0.046) 



-(2.77 + 0.52) 
-(2.73 + 0.39) 
-(2.67 + 0.56) 
-(2.795 + 0.100) 
-(2.774 + 0.086) 

-(2.780 + 0.063) 



Lebrun et al. El 
Brandle et a/.S 
Sugimoto et al^ 
Brandle et al. 
Masuda et al. 
Minamisono et al. 
Minamisono et al. 

Average values 



a large magnetic holding field in which the degeneracy of the Zeeman levels is broken 
by the quadrupole couplings. The polarization is detected through an asymmetry in 
the angular distribution of the electrons in the subsequent [3 decays. 

The first determination^! of a± from the /3-energy dependence of the Pi (cos ^^e) 
term in the decay electron's angular distribution (see Eq. |^) turned out to be fraught 
with experimental difficulties (primarily from backscattering) , since a small effect is to 
be isolated from a large leading term {i.e. aW <^ 1). More reliable data have come 



from later experiments that concentrated on the P2(cos^e) term in Eq. 55. A summary 
of experimental results is given in Table ^. The data are all in accord and the average 
value for the difference in the mirror correlation coefficients is 

Aa = a_ - a+ = 2.737 + 0.078 GeV"^ (56) 
This result is in good agreement with the CVC expectation given in Eq. p2|. 



2.3.4. j3-a correlations 

A similar test of CVC can be mounted in the mass-8 mirror transitions *Li(2"'', 1) — >■ 
SBe(2+, 0) + e" + and 8B(2+, 1) ^ 8Be(2+, 0) + e+ + v^- The broad daughter state, 
which becomes aligned with respect to the electron direction in the process, subse- 
quently breaks up into a pair of a particles. The degree of alignment is measured from 
the cos^ 9 distribution of a-particles, where 9 is the angle between the electron and 
a-particle directions. The (3-a directional correlation in the A=8 nuclei and the elec- 
tron angular distribution from aligned A=12 nuclei are therefore essentially equivalent 
phenomenal with the role of the initial and final states inverted. 

The CVC prediction for this case requires the analogous Ml 7-transition to be 
measured, but this poses more problems than in the mass- 12 case. Firstly, the 2+ 2+ 



sequence implies that the electromagnetic transition is a mix of Ml and E2 multipoles, 
so a mixing ratio has to be measured to determine the Ml fraction. Secondly, the 
analogue state of the decay parent in ^Be is a nearly degenerate doublet (16.6 and 16.9 
MeV) of which both members are roughly equal mixes of T=0 and T=l eigenstates. De 
BraeckeleeJli recently reported on a determination of the Ml and E2 widths needed for 
the CVC test in which the ^}ie{a, 7) cross section was measured as a function of energy 
across the 16.6 and 16.9 MeV levels and photon angular distributions were obtained 
for both resonances. The analysis leads to an isovector Ml width of 2.80 ± 0.18 eV, in 
serious disagreement with previous determinations by Bowles and Garvey0'0. 

The (3-a correlations in the two mirror decays have been measured by Tribble 
and GarveyEll and the results are in agreement with the CVC prediction based on the 
earlier Ml width, although there are large uncertainties. A subsequent experiment^ 
detecting both alpha particles enabled the correlation to be determined as a function 
of the final-state energy. The results show a deviation from the CVC prediction, which 
is dependent on the final-state energy. A repeat experimentEi is currently underway, 
in which it is proposed to measure the triple (5-v-a angular correlation. The neutrino 
direction will be inferred from the recoil of the second a particle. Furthermore, by an 
appropriate choice of geometry, it may be possible to separate the weak-rnagnetism and 
the second-class-current contributions to the alignment. The experimented is expected 
to improve present limits on second-class currents by a factor of five. 



2.3.5. /9-7 correlations 

Another similar test of CVC is based on the mass-20 mirror transitions 20f(2+, 1) 
^ 20Ne(2+, 0) + e" + 17, and 20Na(2+, 1) ^ ^0^6(2+, 0) + e+ + z/,. In this case, the 
daughter nucleus de-excites by 7-ray emission rather than a emission as in the mass-8 
case. The /3-7 angular correlation measurement is equivalent to the l3-a correlation 
experiment and the analysis therefore follows analogously. The result confirm 
the CVC hypothesis, but only to an accuracy of ±25%. 

3. Axial- vector Interaction in the Nucleon 

3.1. Partially Conserved Axial-vector Current Hypothesis 

We start with a small aside and consider the decay of the pion: tt^ i^T^i-, where 
i represents either the first generation, e~z7e, or second generation, \i~V^, leptons. Since 
the pion is a bound state of d and u quarks, its decay is a straightforward example of 
a semi-leptonic charged-current weak interaction. The T-matrix for the decay process 
is given by Eq. ^, and is easily computed for free quarks. However, for pion decay, 
the quarks are bound within the pion and, as mentioned in sect. |1.3| , the required 
quark density functions are not reliably computed. Thus, a more phenomenological 
approach is followed in which we parameterize the matrix element of the quark hadronic 
current, J^, evaluated between the pion and the vacuum, (0| J^|7r~). Since the pion is 



spinless, the only available four-vector is the momentum transfer between the hadron 
and leptons, k^. Thus, the current is parameterized as 



j;: = tUk.TT^ (57) 

where vr" is the field operator for the pion, with the superscript a an isospin index, and 
/vr is some real numerical constant, which parameterizes our inability to evaluate the 
hadron matrix element. In principle, /^r should be calculable in QCD. Note that since 
the pion is a pseudoscalar particle and k^ a polar vector, only the axial portion of J^, 
narnelvA^, is operative in Eq.|5^. The lifetime for the decay is evaluated in a standard 
wayl^3S give: 



T 



--^'i^- ^tvt) = ^yLm\mA\-'^\ . (58) 



The decay is predominately to muons. Thus, with the measured pion lifetime, corrected 
for other small branches and adjusted by a few percent for radiative corrections, a value 
for /j^^is determined from Eq. ^ with the known values of mfj,,m.,r,Gp and Vud- The 



results is 



= 0.663m^ = 92.42 ± 0.26 MeV. (59) 

The constant, is called the pion-decay constant. 

The axial- vector current has no electromagnetic analogue. Furthermore, it is not 
a conserved current since the pion does decay. The axial current and its derivative for 
pion decay are therefore 



Afj, — ifnk^'K"' 

d,Al = ik^Al = -Uky = mlUn\ (60) 

where A;^ = —rn^ for free pions. Its divergence is non-zero because fn ^ and ^ 0. 
However the pion is light on the hadronic mass scale and its mass can be considered 
small. Gell-Mann and LevyS introduced the concept of a partially conserved axial 
current (PCAC) by postulating that conservation is achieved in the "soft-pion limit", 
in which the pion mass tends to zero. Thus, the divergence of the axial current appears 
to be related to the creation operator for a pion. If this principle is also applied to the 
axial current operative in the decay of the neutron, it can be assumed that the current 
is dominated by the pion-pole graph. That is, a pion is created at the proton-neutron 
vertex and then subsequently decays. Schematically, the axial current isH 

^''Note our definition of differs from the Particle Data groupi by \/2 



y4|^ = (n ^ p-R vertex) x (tt propagator) x (tt ^ e vertex) 



^2 _|_ ^2 



■^[^fc;.75]V^^r^ (61) 



where (^j^nn is the pion-nucleon couphng constant. Here F^nn(A;^) is a vertex form factor, 
taken to be a smooth function of k'^ with F.„j^t^{k'^ = —ml) = 1 so that Qt^nn becomes 
the physically- measured coupling constant. Thus, the derivative of this current in the 
soft-pion limit (neglecting the derivative of the vertex form factor) is 

lim = -2 2g^^^F^^^{e)Ui^[-f5]^ \r\ (62) 

Another possible pole graph involves the Ai-meson, a meson of the same spin- 
parity quantum numbers as the axial current. Then the axial current is 



A^^ = [n vertex) x (Ai propagator) x (Ai e i/e vertex) 



"A 




2 



7/^75 + 



V^K' (63) 



where ji = ka'ja- The AiA^A^ coupling constant is extracted from a phenomenological 
chirally-invariant Lagrangianll in which the Ai-meson is treated as the chiral partner of 
the p-meson, and its mass related to that of the p-meson via the Weinberg relationil: 
rriA = V2mp. In particular, we quote from the work of Ivanov and TruhlikP where 
QAiNN = QaQp with Qp the coupling constant for the decay of the rho- meson to two 
pions. Similarly, in the same Lagrangian, the decay of the Ai-meson to leptons is given 
in terms of the same coupling constant, Qp, which then cancels out in forming the axial 
current in Eq. ^ This is the feature of chiral Lagrangians that all vertices involving p 
and Ai mesons are expressed in terms of a common coupling constant, here Qp. Again, a 
smooth vertex function has been introduced with Fp^^^^k"^ = —ml) = 1. The derivative 
of this current is 



dA^ANN^k ) 



mi 



k"^ + m\ 



/^75 + 



k^ /i^75 



mi 



-(7AFANN(r)V'[^75]V'^^" 

2^^Ai"ANN(/c')M^[75]7/;lr^ 



(64) 



with M the nucleon mass. Then, if the pion- and y4i-meson pole graphs, Eqs.^ and 64, 
are added together, the derivative of the axial current in the soft-pion limit becomes 



(65) 



and the requirement that this derivative vanish leads to the condition 



gAFANN{k'^)M = gnNNF^NNik'^)f7r- (66) 

If, in addition, it is assumed that the vertex form factors are slowly varying functions 
of fc^, such that they can be replaced by their limiting values of unity, then 



(67) 



This result is known as the Goldberger-Treiman relation. 

The accepted value for (7,rNN for many years was that determined from 7r-nucleon 
scattering data by Bugg, Carter and Carteril of 13.40 ± 0.08 and confirmed by the 
Karlsruhe-Helsinkiii group. Recently, this value has been challenged, notably by the 
Nijmegen groupS, and smaller values are currently recommended. Three recent deter- 
minations are: 13.03±0.03 (ref.@), 13.14±0.07 (ref.0), and 13.12±0.13 (ref.@). Since 
this issue is still under discussion we will use the last value as a conservative estimate 
for what follows. With this choice for g^NN, the pion-decay constant from Eq.^, and the 
average of the proton and neutron masses, M = 938.919 MeV, the Goldberger-Treiman 
relation predicts for the axial- vector constant | fifA | = 1.29 ±0.01. This is reasonably close 
to the value 1.26 obtained from neutron beta decay (see next section). The difference 
of only 2% is remarkable testimony to the accuracy of the PCAC hypothesis. 

Finally, the sum of Eqs. |61] and yield for the total axial current 



«fi'A-FANN(fc )^ 



^5'a-^ann(^ 



k"^ + ml 



fc2 + ml 



7m75 + 



{mj - ml) 
{k^ + ml){k^ + ml) 
2Mi{l- ml/ml) 



^2 _|_ ^2 



2 ' ' 



(68) 



which can be compared with the phenomenological expression of Eq.[T^ to obtain ex- 
pressions for the form factors ghik"^) and gp{k'^). The results are 



g^ik"^) = gAFAN^ik'^ 



ml 



k'^ + ml 



and 



(69) 



gAk' 



2MgA{e] 

f^2 _j_ ^2 



1 - 



m: 



m 



2 ' 



(70) 



or 



9Ak') = ll^F^^^ik'). (71) 

Since these derivations are based on meson-pole graphs, the form factors deduced 
are of monopole form, with the range given by the mass of the appropriate meson. 
For 5'a(^^)! use of a phenomenological chiral Lagrangian has fixed = \/2mp = 
1086 MeV, a value not too far from the physical Ai meson mas^l of 1230 ± 40 MeV. 
Experimental information has also been obtained from neutrino-nucleon scattering^ 
at the Brookhaven AGS, where the antineutrino quasielastic reaction V^^p — > /i+n was 
studied in the momentum range up to fc^ ~ 1.0 GeV^. A value of an axial- vector mass, 
m^, in a dipole parameterization of the form factor was determined from the shape 
of the momentum distribution to be 1.09 ± 0.04 GeV. This agrees with the monopole 
value, TTT-A, obtained from the phenomenological chiral Lagrangian and suggests the 
vertex form factor Fann(^^) is likewise a monopole function of similar range. 

In Eqs. [7D] and |7T] we gave two expressions for the pseudoscalar coupling constant. 
The first uses the Goldberger-Treiman relation and includes a factor (1 — m^/m^), 
which originates in the propagator for the Ai-meson. The second expression assumes 
that the pseudoscalar coupling constant is entirely due to the pion-pole graph, and 
is the expression we will use here. In semi-leptonic decays, the pseudoscalar coupling 
constant is multiplied by the mass of the lepton, m^, in the expression for the transition 
rate. Thus, for processes involving first generation leptons, such as nuclear beta decay, 
the pseudoscalar term gives a negligible contribution: see sect.^]^. In contrast, for 
processes such as muon capture, involving second-generation leptons, the contribution 
is significant and characterized by the dimensionless quantity 



5'a(0) 5'a(A;^ + ml 

in which ^^^n = 13.12±0.13, = 92.42±0.26 MeV, gj, = 1.2599±0.0025 from neutron 
decay and ^^^^(fc^) ~ 1. A correction to this expression has recently been given by 
Bernard, Kaiser and MeissneiS. It is based on the chiral Ward identity of QCD and 
exploits heavy- meson chiral perturbation theory to find an expression to order fc^. The 
result is 



' 9a{<^) 9A{k' + ml) 3 ^ ^ ^ 



where is the mean square axial radius of the nucleon obtained from the axial-vector 



form factor: 5'a(A; ) = fi'A(0)(l — ^k^rl). In dipole fit^, rj^ = vl2/mo = 0.63 ± 
0.02 fm. corresponding to = 1.09 ± 0.04 GeV. The prediction for the dimensionless 
pseudoscalar quantity then becomes 

ma^^^ = (7.04 ± 0.07) - (0.33 ± 0.02) = 6.71 ± 0.08. (74) 
^a(O) 



In both Eqs. 72 and |7^, the numerical result is obtained for a value of fc^ appropriate 
for muon capture (A;^ = —m'j^ + 2m^Q ~ 0.88m^) where Q is the reaction Q-value. A 
measurement of Qp in muon capture by hydrogen, in addition to providing a test of 
PCAC and the dominance of the pion pole, could, if sufficient accuracy were achieved, 
confirm or deny the correction term predictedii from heavy-meson chiral perturbation 
theory. 

3.2. Neutron Decay 

The decay of the neutron is the simplest (3-decaj mode that involves both the 
vector and the axial-vector interactions. In principle, its study should provide us with 
the most direct view of g^^ and g^, uncluttered by complex nuclear wave functions or by 
renormalization effects of the nuclear medium. In practice, however, the experimental 
uncertainties have, until recently, been large. The value of gv, as determined from 
neutron decay, is a factor of five less precise than the value obtained from the 0^ 0^ 
nuclear /9-decays (see sect.|2.2|). Ultimately, though, its potential is greater since the 
nuclear result has reached a plateau established by theoretical uncertainties that are 
not likely to be soon reduced. 

In the same notation used to describe the superallowed nuclear decays (see 
Eq. the neutron decay rate can be written: 

K 

with 



G'^ = il + Al)Gl 

G'l = (1 + K)Gl. (76) 

Here G^ = GpVud and Ga = GpVudgA with Gp the fundamental weak-interaction cou- 
pling constant determined from muon decay, to which all other coupling constants are 
referenced. For neutron decay, (My)^ = 1 and (Ma)^ = 3. 

Although the calculation of the matrix elements. My and Ma, does not offer any 
theoretical challenge in the case of neutron decay, it is evident that a measurement of 
the neutron's /t-value must be complemented by one additional measurement before 
useful information can be extracted on the individual coupling constants, G'^ and Ga- 
To appreciate what other measurements are possible, we write the differential of the 
probability distribution for neutron (3 decay, n — > pe~z7g, in the following formS: 



d^r oc N{W) 



1 + a-pe ■ Pf + (cr„) ■ A-pe + Bpu 



(77) 



in which N(W) is the usual allowed electron energy spectrum and (cr„) is the av- 
erage neutron polarization vector. The three coefficients, a, A and B, known respec- 



tively as the e-v angular-correlation coefficient, the /?-asymmetry parameter and the 
17- asymmetry parameter, are each related to the ratio of coupling constants: viz. , 

■ ^^-^4^; B-.'StL^, (78) 



1 + 3A2 ' 1 + 3A2 ' 1 + 3A2 ' 

with 

A = ajG, ^ g,. (79) 

All three coefficients have been measured and the ffist two, being more sensitive to A, 
have been combined with the neutron's ft value to determine G'^ and G^. 

3.2.1. /t- value measurements 

As with nuclear decay, the ft value for neutron decay depends upon its Qec value, 
half-life and branching ratio. The Qec value is derived from the evaluated masses^a 
of the neutron and proton, which are known to high precision via a combination of 
Penning-trap mass measurements (see, for example, ref. 0) and measurements of the 
neutron binding energy in deuterium {e.g. ref. S). The result is Qec = 1293.3392(22) 
keV. 

Measurement of the neutron lifetime has defied such high levels of precision. Its 
relatively long, 10.2 min half-life, together with the difficulty in containing a well- 
defined source, have conspired to restrict precision to a half percent or so - more than 
an order-of-magnitude worse than the best nuclear-lifetime measurements - in spite of 
heroic efforts. The measurements (for a recent review, see ref.^) can be separated into 
two distinct categories: those in which the decay products are counted from a defined 
segment of a slow-neutron beam, and those in which cold neutrons are stored and their 
number measured after a well-determined time. These categories have been described^ 
as "counting the dead" and "counting the survivors" , respectively. 

3.2.1.1 Neutron-beam experiments. These experiments depend upon an accurate 
knowledge of the efficiency, e^, with which decay products can be detected, and upon 
the total number of neutrons, A^, being viewed by the decay-product detectors. If 
is the observed decay rate, then 

C, = ^ = (80) 
r r 

Here r is the neutron mean-life (conventionally used in this field instead of the half- 
life: t = rln2), is the neutron density and is the sensitive decay volume. If a 
continuous neutron beam is used, special attention is required to define the boundaries 
of Vd, with a pulsed beam, Vd is defined by the dimensions of the neutron packet, but 
at the expense of a reduced counting rate. 

Results for the five most-precise beam measurements appear in Table 1^. Of the 



three continuous-beam experiments, oneli^ detected the decay electrons, the other 



Table 5. Experimental results'^ for neutron decay 



Method 



Measured values 



Average values 



1 



subsets 



overall 



r(s) n beam 
n trap 

A /3-asyni. 
e-V corr. 



918± 1412^ 



876 ±211 
893.6 ± pi 
903± 13i2: 
877± lOH 
5.4 ±3.3 



104 



loe 



109 



-1.254 ±0.015aLH 
-1.262 ±0.005!ii 
-1.266 ±0.00411^ 
-1.259 ±0.017E! 



891 ±9l 
878 ± 31 



101 



103 



893 ± 2tf 
887.6 ±3.0 
882.6 ±2.7 



m 
ic 



-1.257 ±0.012'liH 
-1.2544 ±0.0036E! 



894.2 ±4.7 



885.9 ± 1.7 



-1.2599 ±0.0026 
-1.259 ±0.017 



887.0 ±2.0 



-1.2599 ± 0.0025 



" Following the practice used for superallowed decay, we retain only those measurements with 
uncertainties that are within a factor of ten of the most precise measurement for each quantity. 
All such measurements not withdrawn by their authors, of which we are aware in early 1995, 
are listed. 

Corrected for weak magnetism and recoil following ref. tl2 



twoE2J'll2a detected protons. The most recent, by Byrne et a/.E23, stored the decay pro- 
tons in a variable-length Penning trap, the magnetic axis of which coincided with the 
neutron-beam axis. To suppress background, the protons were released periodically in a 
brief burst and accelerated into a solid-state detector. The two remaining beam exper- 
iments used pulsed cold-neutron beams: the firstli23 observed decay electrons with an 
in-line superconducting spectrometer, PERKEO; the otherE^l used a helium-filled time- 
projection chamber to observe the decay electrons and measure the neutron density in 
the same apparatus. 

3.2.1.2 Stored-neutron experiments. Neutrons of sufficiently low energy (< 10^''' 
eV) can be successfully confined by material walls, from which they are coherently 
scattered, by gravity or by a magnetic- field gradient acting on their magnetic moment. 
In experiments of this type, a reproducible but not necessarily well known number of 
neutrons is stored and then the number that remains after a storage time, t, is counted. 
The result is related to the neutron mean-life by 



N{t) = N{0) exp 



-t[- + 



1 



^1) 



in which rioss reflects the various neutron-loss processes except (3 decay. The measure- 
ment is then repeated for many different storage times and under different conditions 



that affect the loss processes. Because the loss effects are deemed to be reasonably well 
understood, these results can be extrapolated to give a zero-loss limit, thus yielding 
the neutron mean-life. 

Results for six measurements with neutron traps are listed in Table |: five of the 
experiments used ultra-cold neutrons in gravitational traps with diverse wall materi- 
als; one employed a magnetic storage ring. The earliest trapEsl had aluminum walls, 
adjustable in temperature from 80 to 750 K, and movable plates to vary the neutron 
collision frequency. This same apparatus was later usedllSi with a thin layer of heavy- 
water ice on the walls to reduce losses substantially. A wall coating of hydrogen-free 
diffusion-pump oil was used by Mampe et al^^ to minimize loss in their variable- 
volume chamber, while Nesvizhevskii et a/.Esl chose frozen oxygen on cooled beryllium 
walls. The most recent refinement^ combined oil-coated walls with on-line detection 
of neutrons inelastically scattered by the vessel wall, one of the most important contri- 



butions to neutron loss. Even the magnetic storage ring of Paul et alx^ suffered losses 
because some injected neutrons were captured in unstable orbits; it is believed that 
those neutrons disappeared rapidly and that decay at later times was entirely due to 
j3 decay. 

The lifetime measurements within each category are internally consistent but 
when they are considered as a single body of data the agreement is less satisfactory 
though acceptable: the confidence level is ~ 25%. It has been suggestedlii^ that the 
difference between the average lifetimes for beam and trap measurements could reflect 
the decay n — *• HVe, which would be unobserved in the former. If one takes this 
approach, the averages quoted in Table ^ yield a branching- ratio limit for the hydrogen- 
atom branch of < 2% at the 95% confidence limit. CalculationsEli suggest that the 
branch is actually much smaller (4 x 10~^%), and in deriving the /t- value for the 
neutron /3-decay we have neglected it. 

The statistical rate function with a radiative correction for neutron j3 decay has 
been calculated by Wilkinsonllii. Correcting his result for the current best Qec value, 
and combining it with the overall average mean-life quoted in Table |^, we obtain 

ft{l + 5b) = 1054.4 ± 2.4 s. (82) 
3.2.2. Determination of A 

To date, the most effective measure of A has proven to be the /3 asymmetry from the 
decay of polarized neutrons. Such experiments are designed with a (3 detector on the 
axis of neutron polarization; then spectra are recorded for both polarization directions 
to yield the asymmetry (see Eg. [T^): 

in which iV^ is the electron count rate for each direction of neutron polarization, P is 



the magnitude of the polarization, 6 is the angle between the electron momentum and 
the polarization direction, and v is the electron velocity. In practice, A must also be 
corrected for the effects (~ 1%) of nuclear recoil and weak magnetismEii before it can 
be used to extract a value for A via Eq. [T^ . 

The results of five such experiments are listed in Table ^. Three — li^,Li| rgq^jj^g^j 
electrons to be observed in coincidence with protons in order to reduce background 
eveiits, while the other two obviated the need for proton detection by relying in one 



caseliij on the large volume and high efficiency of the in-line PERKEO spectrometer 



and in the other&i^ on a time-projection chamber operated in coincidence with plastic 
scintillators to determine the electrons' track and energy, respectively. Since the mid 
1980's, the availability of high-flux beams of polarized cold neutrons has helped improve 
the precision of such measurements. 

Only one experiment to determine the e-T7 angular-correlation coefficient, a, has 
achieved sufficient precision to contribute to the determination of A (see Table |^). 
Stratowa et al.^^ measured the energy spectrum of protons recoiling down a tangential 
reactor beam tube from neutrons decaying in the centre of the reactor. The shape of 
this spectrum characterizes the correlation between the unobserved decay products: 
the electron and the antineutrino. 

The five measurements of A based on the observed /^-asymmetry are consistent 
with one another, and with the measurement from the e-V correlation. The resultant 
world average is given in Table ^. 



3.2.3. Implications for the Standard Model 

We are now in a position to extract experimental values for the vector and axial- 
vector coupling constants from neutron decay, compare the former with the result 
from O"*" — > O"*" nuclear beta decays, and examine the implications for the unitarity 
test of the CKM matrix (see sect. |2.2.^ . The current experimental status from neutron 
decay can be summarized in two expressions, 



{G'^ + ?>G'l)/{hcf 



-1.2599 ±0.0025 
(7.701 ±0.018) X 10"^° 



GeV" 



where the latter follows from combining Eqs.|7^ and g2[ The solution of these equations, 
shown graphically in figure 2, gives the following results for the effective coupling 
constants: 



G'J{hcf = -(1.4566 ± 0.0018) X 10"^ GeV^^ (85) 
G'^/{hcY = (1.1561 ±0.0023) X 10"^ GeV"2 [neutron]. (86) 




Fig. 2. Allowed regions of and G'y from neutron decay, nuclear superallowed decays and 
CKM unitarity 



In sect. |2 .21 we described the study of 0^ — > 0^ /3-decay transitions in nine nuclei. 
The results, shown in Table 2, led to an average J^t value which, with the help of 
Eqs. |2^ and B^, becomes: 



G'^/ihcf = (1.1494 ± 0.0006) x 10"^ GeV 



[0^ 



0' 



(87) 



where the quoted error bar includes provision for a "systematic" uncertainty in cal- 
culating the charge corrections (see sects. 2.2.3.1 and p.2.4|). The neutron result is 



considerably less precise but, even so, it is in serious disagreement with nuclear decays. 
The discrepancy is not understood. The analysis of the nuclear results suffers from un- 
certainties in the effects of charge dependence and nuclear structure; the neutron work 
is faced with more formidable experimental challenges. There is no reason to suppose, 
at this time, that either result is in doubt beyond its quoted uncertainty. 

Following the approach taken in sect, p.2.^ the neutron result for the vector 
coupling constant can be used to test the unitarity of the CKM matrix, yielding 



I Vud I '+ I Vus I '+ I Kfe I ' = 1.0077 ± 0.0040, (88) 

which differs from unitarity at a confidence level of ~ 99%, but in the opposite sense 
to the result from nuclear (3 decay. In the absence of any viable explanation within the 
Standard Model for the discrepancy between the neutron and nuclear results, nor any 



reason to prefer one over the other, we can do no better than take the weighted average 
of both results for Gy, suitably inflating the resultant uncertainty. This leads to 



G'^/{hcf = (1.1498 ± 0.0015) x 10"^ GeV"^ [average] 

iKdl = 0.9741 ±0.0017 [average] (89) 



and 



I Kdl^+lKs I ^+1 Kfel^ = 0.9978 ±0.0037 [average], (90) 

which agrees with unitarity. 

3.3. Muon Capture in Hydrogen 

The induced pseudoscalar couphng constant, g^, of Eq.|TD|, is the most poorly 
determined of all the nucleon weak-interaction coupling constants. This is because in 
semi-leptonic weak processes involving first-generation leptons, such as beta decay or 
electron capture, the pseudoscalar term gives negligible contribution to the transition 
rate. Consider the T-matrix involving just the pseudoscalar term: 



G — — 

Tfi = -^Kd#p(i5p^M75)V'n^V'e7M(l + 75)^^77^ 



G 

= -;^Kd#p(^^ ^75)^n^e(l + 75)^/^77^, (91) 

where k = pe+Vu- Here pseudoscalar coupling is transformed to pseudovector coupling 
in line two, a step valid only for on-shell nucleons. The factor /c^ is moved into the lepton 
matrix element in line three, where application of the Dirac equation, = imeipe^ 

shows immediately that the pseudoscalar term is of size gpme/[2M) and is very small. 
However, for second-generation leptons participating in such processes as muon capture, 
the factor becomes gpm^/ (2M) and is much larger. 

Consequently, experiments to determine g^ are based on muon capture in hydro- 
gen. The two main processes that contribute to the disappearance of a negative muon 
stopped in pure hydrogen are muon decay, /i~ — > e~ ± I7e ± i/^, and nuclear capture, 
/i~ +p n + u^. In the latter process, a 5.2 MeV neutron is emitted. Because the weak 
interaction is strongly spin dependent, an accurate specification of the initial muon- 
proton state is vital to the interpretation of the experimental results. The history of a 
negative muon, stopping in a hydrogen target, is as follows. At low velocity, the muon 



is captured by a proton and cascades to the IS" state in less than 10~^s; this neutral 
(hp) atomic system is formed in a statistical mixture of triplet and singlet states. In 
general, the triplet state converts to the singlet state with a rate depending upon the 
hydrogen density. Even in gaseous hydrogen, for pressures above five atmospheres, the 
conversion rate is large compared to the natural decay rate, so nuclear capture occurs 
predominantly from the singlet state. 

In liquid hydrogen, the singlet (/ip) system further converts rapidly into a {pfip)~^ 
molecular ion if the muon has not decayed first. Consequently, it is from the stable 
molecular ion that the muon is most likely to be captured or to decay. Most experiments 
are able to measure the time of each capture and thus the configuration from which 
capture takes place is known. Usually only capture events that have originated in a 
[pfip) system are included in the analysis. The {pfip)~^ molecular ion is almost invariably 
formed in the ortho state for which the capture rate, Aom, is given by 27o(|As + \At), 
where A refers the the capture in the triplet and singlet atoms, and 270 is a wave 
function overlap factor of the muon with the proton, calculated to be 27^ = 1.01 ±0.01. 
There is a small probability that the ortho state will decay to the lower-energy para- 
molecular state before the capture occurs. The capture rate from the para state is given 
by Apm = 2jp{jAs + |Aj), with 27^ = 1.15 ± 0.01, and is smaller than the ortho rate 
because A^ ^ Aj. Nevertheless, a correction to the measured capture rate has to be 
applied for this ortho-to-para transition and this correction significantly influences the 
extraction of the gp coupli ng constant. 

Early experimentJi^^'SEHi'EHl obtained the capture rate by observing the out- 
going 5.2 MeV neutrons. The principal limitation was the calibration of the absolute 
efficiency of the neutron detectors. A more recent Saclay experimentli^ obtained the 
capture rate by comparing the lifetime r^- of negative muons stopped in liquid hy- 
drogen with the lifetime of positive muons. Both r^- and were measured by 
detection of the delayed muon-decay electrons. The capture rate was then obtained 
from the difference in these two lifetimes. In all cases, the hydrogen in the target had 
to be highly purified and isotopically pure. If small amounts of deuterium or other 
atoms are present, then an irreversible transfer to the deuteron (or heavier atom) is 
likely. Use of pure hydrogen, therefore, is very important to minimize this loss. The 



Columbia and Saclay experimentalH-Jf H I used liquid hydrogen targets; the capture 



rates were corrected for the ortho-to-para transition, which was explicitlv measured in 



a subsequent Saclay experimentllH3. The Dubna and CERN experimentsllH3'EH3 used gas 
targets. 

The resulting values of the dimensionless quantity m^gp/o^ from these five exper- 
iments are: 4.8±6.3 (ref.El), 8.7±3.7 (ref.El), 8.2±3.1 (ref.ti), 6.3±4.7 (ref.lll) and 
5.6 ± 2.6 (refs.li^lllHi). The average is 6.9 ± 1.5 in excellent agreement with the PCAC 
value of 6.7 from Eq.|7^, evaluated with k"^ = O.SSm^, which is appropriate for hydro- 
gen. Were it not for the ortho-to-para correction, the value from the Saclay experiment 
would have been 10.6 ± 1.6 rather than the quoted 5.6 ± 2.6! A new measurement of 
the ortho-to-para transition rate with improved accuracy would be very desirable. 



Despite the good agreement between the ordinary muon capture (OMC) experi- 
ments and the PCAC prediction, it should be noted that the OMC rate is insensitive 
to the exact value of Qp since an individual measurement, based on a 4% OMC rate 
determination at Saclay^, gives a 42% uncertainty for gp. This is because, in liquid 
hydrogen, 80% of the capture takes place from the ortho-molecular state, which is pre- 
dominantly the atomic singlet state. In contrast, in radiative muon capture (RMC), in 
which OMC is accompanied by a bremsstrahlung photon, the photon carries away one 
unit of angular momentum and the capture is mainly from the triplet atomic state. 
The triplet capture rate is 15 times more sensitive than the singlet rate to the value 



of Qp. Beder and FearingllHj have estimated that an 8% measurement of the RMC rate 
could provide a 10% value for gp. 

The experimental problems in measuring RMC in hydrogen are formidable. The 
partial RMC branching ratio (for > 57 MeV) is only 1.6 x 10~^ and the bremsstrah- 
lung from radiative muon decay or from the Michel electrons is nearly five orders of 
magnitude larger. Consequently, the apparatus requires good energy resolution and an 
accurate knowledge of the high-energy tail of the response function in order to guarantee 
that the RMC signal region is not contaminated by events from the bremsstrahlung 
region. The low branching ratio also requires an extremely high purity level in the 
hydrogen liquid itself; otherwise the fi~ will be lost from the proton before the capture 
can take place. Impurity levels of < 10~^ for all elements other than hydrogen are 
necessary since the capture probability increases as Z^. An experimentHZl presently 
underway at TRIUMF hopes to achieve an 8% measurement of the RMC capture rate. 

4. Axial-vector Interaction in Nuclei 

4-1- Quenching of the axial-vector current in nuclei 

Since the vector current is a conserved current, its normalization is unchanged 
from one environment to another: g^iO) is unity at the quark level, at the single- 
nucleon level and in finite nuclei. In contrast to this, the axial current is not conserved 
and its normalization therefore becomes environment dependent: (y'A(O) is unity at the 
quark level, is greater than unity (of order 1.26) at the single-nucleon level and reduced 
relative to 1.26 in finite nuclei. The degree of reduction, however, is nucleus dependent. 
Since pions cause the failure of the axial current to be conserved, as evidenced in 
the PCAC hypothesis (see Eq. ^Up, and since nucleons in nuclei are interacting by the 
exchange of pions, it is not surprising that the axial current in nuclei is modified. 
One mechanism frequently discussed in this context involves the A(1232, s = 3/2, t = 
3/2)-isobar, which can be viewed as an internal spin-isospin-flip excitation of a single 
nucleon. A nucleon, for example, in emitting a pion can be raised to its A-excited 
state. The A is subsequently de-excited back to the nucleon through an interaction 
with the axial current, while the pion is absorbed by a second nucleon. This and 
related processes are calculable in perturbation theory but the results are strongly 
model dependent, relying on assumptions made concerning the isobar's interactions. It 



is probably worth repeating here the cautionary remarks of Lipkin and Leet=^: "There 



has been considerable debate of the possibility that isobars may be present in nuclei. 
However, it is not clear exactly what this means. The isobar is a highly unstable 
particle with a natural width much larger than the spacing between nuclear levels. It 
is, therefore, unclear whether an isobar present in a nucleus can be considered as an 
elementary fermion in some approximation, as a three-quark composite, or as a pion- 
nucleon resonance. Furthermore, an isobar present in a normal nuclear ground state 
is very far off shell, and it is not clear how to extrapolate the properties of such a 
comphcated object far off shell. All this confuses the issue and makes it difficult to test 
experimentally for the presence of such isobars." 

The idea that the axial-vector coupling constant might be renormalized for nu- 
cleons embedded in a nucleus was first discussed by Ericson et al. 1^23. They pointed out 
that the pion field within the nucleus might be suppressed at small momenta because of 
the strong short-range repulsion between nucleons in the medium. Thus, an axial-vector 
/?-decay in nuclei would, using the arguments of PCAC, be similarly suppressed. This 
was worked out in more detail by RhoEl and Ohta a„d WakamatsuEl, who showed 
that the A-isobar would be the main mediator in effecting such a suppression. Exper- 
imental data, as we will shortly discuss, also give evidence for this suppression, but 
whether or not it is due to the influences of A-isobars is not so clear. There are other 
many-body corrections to be applied to the nuclear-structure calculations, which arise 
from short-range central and tensor correlations between nucleons and which also lead 
to a suppression of axial-vector ratesliH'^. 

One of the primary sources of experimental information on axial- vector transitions 



in nuclei comes from nuclear beta decay. Referring to the axial current given in Eq. |T2 
we ignore second-class currents and note that the pseudoscalar coupling term gives 
a negligible contribution. Thus, only the iipgAlfil^^T""^ term is operable. In nuclear 
physics it is sufficiently accurate to consider a nonrelativistic limit of this current in 
which only leading terms in p/M, where p is the nucleon momentum and M its mass, 
are retained and four- component Dirac spinors, ip, are replaced by two-component 
Pauli spinors, x- The leading term arises when the Lorentz index, fi, is space-like: 

A'^ = -gj^x^fT^T^x, (92) 

where cr is the Pauli spin operator and = ± ir^, with the upper sign for /3~ 
transitions and the lower sign for /3+ transitions. The quantity between x^X is taken 
as a one-body operator, the matrix elements of which are evaluated in the many-body 
nucleus with shell-model wavefunctions. The experimentally determined quantity, the 
/t- value, for a pure axial- vector transition is related to this matrix element via (c/. 
Eq.m 



where 



(93) 



B{GT) = {M^f 
(M^) = -^{T,T,i 1 T 1 I TfT,f){JfTf III CTT III J.Ti). (94) 

Here the spin, isospin and z-projection of isospin are denoted Ji,Tj and T^i for the 
initial state and Jf,Tf and Tzf for the final state. The triple-barred matrix element is 
a reduced matrix element in both spin and isospin spaces in the notation of Brink and 



Satchlerfi^ and J = (2 J + 1)^/^^ The constant, 6146 ± 6 s, in Eq.|9| is given by 2 J't, 
where JFt is the average electromagnetically-corrected ft value for the superallowed 
Fermi O"*" — >■ O"*" beta transitions (see Eq.|39|). 

The operator err, known as the Gamow- Teller (GT) operator, has no dependence 
on orbital angular momentum and hence has no non-zero matrix elements between 
states of different orbital symmetry. In particular, in a shell-model calculation carried 
out in a complete basis of harmonic oscillator states {e.g. a (0p3/2, 0pi/2)-basis for p-shell 
nuclei or a (0^5/2, lsi/2, 0(i3/2)-basis for sci-shell nuclei) the GT-operator acting on any 
member of the basis connects only to other states within the basis and not to any state 
with occupancy in a different major oscillator shell. Complete untruncated shell-model 
calculations within a single major oscillator shell are called Ohu calculations. A working 
definition for the term quenching is as follows: for beta decays populating well-defined, 
isolated states in the daughter nucleus, the ratio of the experimentally measured rate to 
the calculated rate in a full Ohu calculation is called the quenching factor. An average 
quenching factor implies an average over many transitions. The definition is model 
dependent in the sense that the shell-model calculation is dependent on the choice of 
residual interaction, but for an average value it is hoped that this dependence is not 
severe. Equivalently, one can discuss a quenching in terms of a renormalization of 
the axial- vector coupling constant. We define 

9a,oh = 9a + SgA, (95) 

where (/a is the free-nucleon value of 1.26 and 6gj^ the correction to it; thus, the quench- 
ing factor is {gA.,.s/9Af- 

The first systematic attempt to determine an average quenching factor, g, was by 
Wilkinson^ from a study of /^-transitions in the p-shell and the lower S(i-shell nuclei 
(yl=17-21), in which he obtained q = 0.80 ± 0.06. A more comprehensive analysis by 
Brown and Wildenthallil using full Ohu calculations throughout the sd-shell for 189 
data, for which the experimental error on the deduced GT matrix element was 10% 
or less, found an average quenching factor of g = 0.63 ± 0.01; the quoted error here 
is limited to the statistical error on the fit. Significant quenching seems to be evident. 
One configuration has a special role to play here, for example in the ground states of 



T^^However, see some contrary remarks by BrownlH3 in a discussion of the /3-decay of '^^Ca 



A = 15, 17, 39 and 41 nuclei: it is a closed-LS-shell-plus (or minus)-one nucleon. In 
these cases, the full Ohu calculation involves just a single-particle matrix element, so 
microscopic calculations can seek to evaluate corrections to the single-particle matrix 
element and identify the origin of the quenching phenomenon. The TokyoE^'liii and 
Chalk Riveiiill^ groups have attempted this with similar results. 

Generally, the quenching stems either from the inadequacy of the single-particle 
description of the nuclear state or from the inadequacy of the one-body GT operator. 
Calculations of the first effect, frequently called core polarization, estimate through 
perturbation theory the admixtures of 2p-lh and 3p-2h configurations in the basi- 
cally single-particle wave function and evaluate the impact the admixtures have on the 
calculated matrix element. At LS closed shells there is no contribution to the core po- 
larization in first order because the one-body GT operator has no spatial dependence 
and thus cannot excite a Ip-lh state from a closed LS shell. Therefore, calculations 
must be taken to second order in perturbation theory. This leads to a time-consuming 
calculation, as there is no selection rule to limit the intermediate-state summation and 
the convergence is slow. This is particularly true with tensor forces in the residual 
interaction as first stressed by Shimizu et al.^^. This propensity of the tensor force 
to couple strongly to high-lying states has been called 'tensor correlations' and the 
phenomenon leads to a reduction in the GT matrix element. 

The other inadequacy concerns the use of one-body operators. Corrections arise 
because nucleons in nuclei interact through the exchange of mesons, and this exchange 
can be perturbed by the action of the weak axial current. Since this perturbation 
requires at least two nucleons to be involved, corrections of this type lead to two- 
body GT operators. The possible pion-exchange processes can be classified into two 
types: Born graphs, involving only nucleons and pions, and non-Born graphs in which 
the axial current either excites a nucleon to its isobar A state, or converts the pion 
into a heavy meson, such as the p. The important point is that for axial currents the 



pion Born graphs are identically zerotlH. For the non-Born graphs, we can distinguish 
between diagrams involving nucleon excitations - referring to these as isobar graphs - 
and those in which a pion is converted to a p meson - referring to these as pvr or meson- 
exchange-current (MEG) graphs. In general, the MEG graphs give a small contribution 
and are not an important ingredient in the quenching of the GT matrix element. 

There are two further points to consider. The core-polarization calculation cor- 
rects the matrix element of a one-body operator evaluated in the closed-shell-plus-one 
configuration for the presence of 2p-lh and 3p-2h admixtures in the single-particle wave 
function. The perturbation calculation is carried out to second order in the residual 
interaction (or to the fourth power in the meson- nucleon coupling constants). It is 
logical, therefore, that the matrix elements of two-body operators should likewise be 
corrected for 2p-lh and 3p-2h admixtures. Since the two-body operator itself involves 
the meson-nucleon coupling constants to the second power, it is sufficient to estimate 
this correction to first order in the residual interaction. These terms have been called 
'crossing terms' in the work of the Tokyo group^. The second point concerns the one- 



Table 6. Corrections to the ground-state diagonal and spin-flip ofF-diagonal matrix elements, 
5{Ma)-, expressed as a percentage of the single-particle value, {Ma)sp-, for closed-shell-plus-one 
configurations at ^ = 17 and A = 41. Quenching factors, g, are also given. 



6{Ma)/{Ma)sp in % 





"5/2 "5/2 


"5/2 "3/2 


/7/2 /7/2 


/7/2 ~^ J5/2 


Core polarization 


-9.0 


-11.3 


-11.8 


-14.1 


Isobars 


-0.9 


-5.9 


-2.1 


-6.2 


MEC 


-1.2 


0.3 


-1.4 


0.2 


Crossing terms 


3.0 


2.3 


3.4 


2.6 


Relativistic corr 


-2.1 


-1.2 


-2.3 


-1.3 


Total 


-10.2 


-15.8 


-14.2 


-18.8 


Quenching factor 


0.81 


0.71 


0.74 


0.66 


Experiment, q 


0.74 ± 0.01 




0.54 ± 0.01 





body GT operator, which is obtained as a leading term in the nonrelativistic reduction 
of a relativistic axial- vector current. There is a few per cent correction coming from 
the next-order terms in the reduction. 

In summary, then, corrections to lowest-order shell-model estimates of the GT 
matrix element come from: core polarization, isobar currents, MEC currents, crossing 
terms and relativistic corrections. All these ingredients for closed-LS-shel l-pl us-one nu- 



clei have been calculated by Arima et al. and by Towner and Khannat^. We quote 



in Table 6 some results from a review by Townertid, in which more details can be found. 

Note that the calculated quenching factors, although in the expected range of 
0.66 to 0.81, are higher than the experimental values for the diagonal matrix elements 
and, for the Od orbits, are also higher than the average quenching factor of 0.63 ± 0.01 



determined by Brown and Wildenthalt2II. Furthermore, there is more quenching evident 
in spin-flip than in diagonal matrix elements, with approximately two-thirds of the 
quenching coming from core polarization (notably tensor correlations) and one-third 
from isobars. 

4-2. {p,n) and {n,p) reactions 

The reactions {p, n) and (n, p) are comparable to 13~ and (3^ decay, respectively, 
but are not limited by the energetics of a radioactive decay. In studying these reactions, 
it is useful to define a Gamow- Teller strength function by the relation: 

S{u) = E(^a(^ - f))'S{Ef -E.-oo), (96) 
/ 

where u is excitation energy, S is a resolution function, and i and / represent the initial 



and final nuclear eigenstates of energies Ei and Ef. The major peaks in S{uj) are iden- 
tified with a collective state: the Gamow- Teller giant resonance (GTR). Unfortunately, 
beta decay has access to nuclear states in a very limited energy window and misses the 
strongest states in the strength function. Therefore, to map out the complete strength 
function one needs probes that allow the independent variation of both energy transfer 
and momentum transfer to the target. This is possible with hadronic probes such as 
the {p, n) and (n, p) charge-exchange reactions. 

The Gamow- Teller resonance has quantum numbers J'^ = 1^ , L = 0, S = 1,T = 1 
indicating it is populated through spin- and isospin-excitations. In the {p, n) reaction 
this mode is reached via the V^rCp ■ f^j^p ■ Tj component of the interaction between 
the projectile and target nucleons. The labels p and j refer to the projectile and the 
struck nucleon respectively. With respect to the target, the (p, n) probe has a similar 
spin-isospin operator structure as the GT operator of /9-decay. However, the (p, n) 
reaction simulates the /5-decay conditions only if the cross-section is measured at very 
small momentum transfer. In the (p, n) reaction, this condition can be met only for 
0° scattering and at high bombarding energies. As there is no orbital excitation in a 
GT transition, the angular distribution for the (p, n) reaction to the GTR will have a 
characteristic L = shape, with the maximum cross-section at zero degrees. 

The GTR occurs at a higher excitation energy than the isobaric analogue of the 
target ground state (IAS), which carries the quantum numbers J'^ = 0+,L = 0,5' = 
0,T = 1. The {p,n) reaction excites the IAS through the force component VrTp ■ Tj, 
which requires no spin flip [S = 0). While the IAS is the dominant peak in the spectrum 
at low incident energies {Ep < 100 MeV), it is completely swampedEll by the large GT 
cross-section at Ep = 200 MeV. The large difference in the excitation strengths of 
the IAS and GTR as a function of the incident projectile energy is due to the strong 
energy dependence of the V^-term in the projectile-target-nucleon interaction. The force 
strength is greatly reduced with increasing projectile energy, while the force strength 
Vf^r is nearly independent of incident energy. 

The {p, n) reaction is a powerful spectroscopic tool for measuring the GT strength 
functions of nuclei. The zero-degree {p, n) cross-section can be calculated in a distorted- 
wave theory and compared to experiment to extract the spectroscopic information: 
namely B{GT) of Eq. 0. However, there is sufficient uncertainty in the distorted- 
wave theory, especially with regard to its absolute normalization, that it is usual to 
calibrate the theory by measuring the {p, n) cross-sections of specific transitions for 
which the B{GT) values are already known from /3-decaytl3. Once this is done, then 
the [p, n) reaction can be used to determine the GT matrix elements of states that 
cannot be populated by /5-decay. In particular, it is possible to map out the whole 
GT response function in the zero-degree spectrum. There are, however, other states 
with multipolarities L > 0, which also contribute to the 0° (p, n) cross-section. In 
general, these states receive their maximum cross- sect ions at larger scattering angles, 
but their angular distributions extend forward to 0°. These states lead to a background 
contribution that has to be subtracted from the data. A discussion of the theoretical 



techniques for doing so can be found in the papers of Osterfeldlll3'tl3. 

The breakthrough in the use of the {p, n) reaction to study GT resonances came 
in 1980 from e xpe riments at the Indiana University Cyclotron Facihty. These experi- 



ment sti3^i^[!^3lll3ll£y demonstrated the existence of very collective spin-isospin modes 
in nuclei. An excellent review of these and the many subsequent experiments is to be 
found in ref. Ell. The spectra of nuclei with neutron excess are seen to be dominated by 
one prominent peak, which is interpreted as the giant Gamow- Teller resonance. This 
collective mode was already predicted by Ikeda, Fujii and FujitaEIll, who inferred the 
existence of the GTR from the absence of spin-isospin strength at low excitation ener- 
gies. Their argument is based on a sum rule. Let us integrate the strength function of 
Eq. ^ to define a sum: 

S = [ Siuj)dLU = J2{M^it ^ f))'. (97) 
/ 

Then, from the commutator algebra of isospin, [r^,r^] = 2zr^, we obtain for the dif- 
ference in the sum rule for /5~ and (3^ emissions from a given nucleus the model- 
independent resulti^: 

S {(]-)- S{f3+) = 3{N - Z). (98) 

That is, the difference in the sums for the {p,n) and {n,p) reactions on a particular 
target is quite independent of nuclear structure. In heavy nuclei with large neutron 
excess, S{/3~^) ^ 0, because all states into which a proton could be transferred within 
the same major shell are already fully occupied: the GT transition is Pauli blocked. 
Setting S{(3^) = then gives a lower limit for the sum: viz. 

S{f3')>3{N - Z). (99) 

At intermediate energies {Ep ~ 100 — 200 MeV), the selective excitation of the GT giant 
resonance by the (p, n) reaction has been observed in almost the entire periodic table. 
The striking resulttla is that experimentally only (60 ± 10)% of the expected sum- rule 
limit, 3(A^ — Z), can be located in the excitation region where the major GT peaks 
occur. This means that at least 40% of the GT strength is unobserved in the spectra. 
As described in the previous section for beta decay, this missing GT stren gth has been 



interpreted as evidence for isobar admixtures in the nucleonic wave function 



12 



130 



I or, 



alternatively, as indicatingt^ that high-lying 2p-2h states mix with the low-lying Ip-lh 
GT states and shift GT strength above and below the GTR. If this latter mechanism 
is the right one, then the so-called missing GT strength would actually be located in 
the background. The detection of this strength would not be an easy task because the 
spectrum on the high-energy side of the GTR is continuous and structureless, making 



it difficult to extract information. Calculations by Smith and Wambacht^ indicate that 



significant strength from 2p-2h mixing resides in the tails of the resonances, suggesting 



that the quenching due to A-isobars is 10 — 20%. This is a further indication that 
quenching due to A's is small. 

In conclusion we note two recent experiments that have elicited an exchange of 
comments in the literature: 



A groupll^ at TRIUMF has measured angular distributions of the ^^Fe(p, n)^^Co 
and ^^Fe{n, p)^^Mn cross-sections at incident beam energies of 300 MeV to test 
the GT sum rule, Eq. ^ for a case in which the GT strength is strong in both 
channels. The {p, n) cross-section at 0° for L = was summed up to an exci- 
tation energy of i^^^ = 15 MeV. The strength above this excitation energy was 
not included because of increasing uncertainties in removing dipole and higher 
multipole strength at higher energies. The integrated strength below = 15 
MeV is S{(3~) = 7.5 ±1.2 in agreement with earlier resultal^ at 160 MeV of 
S{l3~) = 7.8 ± 1.9. The 0° {n,p) cross-section was summed up to = 10 MeV, 
the same criteria being used to decide the upper limit of the integration. The 



integrated strength is S{f]^) = 3.1 ± 0.6. A more recent experiments^ at a lower 
energy of 97 MeV finds 3.5 ±0.3 ±0.4, the first uncertainty being due to statistics 
and multipole decomposition, the second due to uncertainties in normalization. 
The sum rule difference in the TRIUMF data thus yields 

S{(3-) - S{(3+) = (7.5 ± 1.2) - (3.1 ± 0.6) = 4.4 ± 1.3, (100) 

to be compared with the model- independent value of 3(A^ ~ = 6 in this 
case. The result is consistent with the quenching factor of 0.6 seen throughout 
the periodic table in {p, n) reactions where the {n, p) strength is strongly Pauli 



blocked. However, Goodman et al. 113 have commented that this work does not 
test the model-independent sum rule because both measured strength functions 
are incomplete, there being additional strength undetected either because it lies 
outside the measured energy domain or because it is hidden in the background. 



A recent beta-decay experiment by Garcia et al. has extended the measured 
GT strength function in the (3^ decay of ^''^Ca up to 8 MeV excitation in the 
daughter nucleus. Delayed protons were detected in coincidence with gamma rays 
to identify proton decays of ^''K daughter levels to excited states in ^^Ar. Some 
43 previously undetected proton groups were observed yielding an integrated 
B{GT) twice as large as that observed previously. In terms of the sum-rule limit, 
the experiment observed a summed GT strength of 1.95 compared to the sum- 
rule requirement of at least 9 units. Because ^''Ca is the isospin mirror of ^^Cl, it 
is possible to compare the strength seen in this decay with that measured in the 



^^Cl(p, n)^'^Ar reaction studied by Rapaport et a/. tEl. Garcia et al. made such a 
comparison and concluded that the decay contains 50% more GT strength in 
this energy region than the (p, n) experiment. 



This conclusion led to a flurry of correspondencal^t^. Two issues were raised. 
First, in the analysis of the {p, n) data, insufficient Gamow- Teller strength was 
assigned near the isobaric analogue state, which in the reaction is populated by 
both Fermi and GT operators. Second, the beta-decay data for the low-lying 
states have to take account of both delayed-gamma and delayed-proton emis- 
sionli 
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The current situation is that the ratio of B{GT) values determined from beta- 
decay and (p, n) experiments show a state-by-state scatter with discrepancies 
in the B{GT) values that can vary by as much as a factor of two. However, 
these states are below 8 MeV excitation and their B{GT) values are very small. 
Discrepancies of this size, therefore, are not necessarily surprising because, for 
weak transitions with B{GT) ^ 1, one can expect inadequacies in the reaction 
model used to extract B{GT) values from forward-angle charge-exchange cross- 
sections. The importance of the A = 37 example is that it allows a comparison of 
GT strength inferred from {p, n) experiments and beta-decay experiments over a 
significantly larger excitation energy range and for significantly more transitions 
than any other system studied so far. Further experiments are underway. 

Despite the large number of GT transitions identified in this example, the (5^ 
experiment found only 2/9 of the sum-rule limit. BrownEH points out that if 
experiment and theory are compared over an energy region that only exhausts 
a small fraction of the sum rule, the results are very sensitive to the choice of 
residual interaction used in the shell-model calculations. Some older interactions 



reproduced the data better than the preferred Wildenthal interactions. This, 
claims Brown, indicates a need for further improvement in the sd-shell Hamilto- 
nian. 



4-3. Radiative Muon Capture in Nuclei 

It was explained in sect. ^]3| that muon-capture data provide the best opportunity 
to measure the induced pseudoscalar coupling constant. Data on ordinary muon capture 
(OMC) by the proton yield a value consistent with the PGAC hypothesis, but the 
accuracy is not sufficient to provide a demanding test of the prediction. Radiative 
muon capture (RMC) potentially provides greater sensitivity, but the low total capture 
probability in hydrogen (~ 10~^) combined with the low expected branching ratio 
RMG/OMC (~ 10^^) has so far precluded a measurement of the elementary radiative 
capture process, although results from an experimentli^ at TRIUMF are anticipated 
soon. 

RMC measurements are easier for heavier nuclei, in which the muon capture 
probability is greatly enhanced: roughly by a factor of Z^, where Z is the charge 
number of the capturing nucleus. Exclusive nuclear RMC to a specific final nuclear 
state has not yet been attempted and would pose formidable experimental challenges. 
All measurements of RMC to date have been inclusive ones, in which the final nuclear 



states are not experimentally resolved. These suffer the disadvantage that the nuclear 
response has to be calculated from some theory and this introduces considerable model 
dependence. 

Most nuclear RMC calculations use the impulse approximation, in which the 
amplitude for the elementary process is summed incoherently over all the nucleons in 
the nucleus. Even within the impulse approximation, the calculation is non-trivial, since 
the RMC amplitudes must be summed over all final states. The earliest calculations 
adopted the closure approximation l^'SE^^ for this final-state summation, in which a 



single average excitation energy with a corresponding maximum photon energy, fcmax, 
is assigned to the final nucleus. The RMC photon spectrum then takes a polynomial 
form, 

(l-2x + 2x2)x(l-x)^ (101) 



where x = E^,/krr,«^ and is the photon energy. The use of the closure approximation 
has been criticizedt^a and, indeed, such calculations typically overestimate the RMC 
branching ratio by large factors^. Nevertheless, the polynomial form given in Eq.|lUl 



usually reproduces the data fairly well provided fcmax is treated as a variable parameter. 
Christillin, Rosa-Clot and ServadioS were able to avoid the closure approximation for 
medium-heavy nuclei by integrating over the nuclear response obtained in a nonrela- 
tivistic Fermi-gas model. The only free parameter was the effective nucleon mass, M*, 
which was fixed at M* = 0.5M from a global fit to experimental OMC rates. The RMC 
rate was then calculated in a consistent manner. A recent calculation by Fearing and 



Welshlllll cast some doubt on the validity of the nonrelativistic Fermi-gas model cal- 
culations. They adopt a fully relativistic Fermi-gas model along with the local-density 
approximation and realistic nuclear-density distributions and a number of other refine- 
ments; they find values for the OMC and RMC rates that are different by factors of 
two or more. This issue remains unresolved. 

In lighter nuclei, a more microscopic account of the response function can be at- 



tempted. Christillinll^ uses a phenomenological nuclear excitation spectrum for ^°Ca 
consisting of Lorentzian giant dipole and quadrupole components. The widths and rel- 
ative strengths are adjusted to reproduce existing photoabsorption and OMC data. A 
similar approach is used by Christillin and Gmitrolili for ^^O except explicit account 
is taken of low-lying states in the ^^N final nucleus in addition to the giant dipole and 
quadrupole terms. Lastly, Roig and Navarroli^ have invoked SU(4) symmetry and ap- 
plied sum-rule techniques to calculate RMC on ^^C, ^^O and '^°Ca. The branching ratios 
predicted are substantially lower than those provided by either the phenomenological 
or microscopic models. 

In Table 7 we tabulate the branchi ng r atio and Qp values extracted from RMC 
measurements at PS10 and at TRIUMF0Jil3 with various theories for the nuclear 



response functions. It is evident that the observed branching ratio is falling for heavier 
nuclei. It also appears that Qp is being quenched but this result is heavily dependent 



Table 7. Branching ratios and pseudoscalar coupling constants extracted from RMC measure- 
ments at PSI and at TRIUMF with various theoretical calculations of the nuclear response 
function. 



Branching 
Target Ratio 



12 



16 



C 



o 



27A1 

40Ca 



Mo 
Sn 

165jjo 

Pb 

209Bi 



m 



CRS80'* 
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on the theory of the nuclear response, which is unsatisfactory at the present time. 

In light nuclei, there are a few OMC experiments to explicit final nuclear states 
that lead to a determination of the pseudoscalar coupling constant. To isolate the 
pseudoscalar term from a term involving the regular axial- vector current requires either 
two complementary measurements on the same transition or one measurement and a 
calculation. Further, with the pseudoscalar term identified, a calculation of the matrix 
element is required to determine the coupling constant. We cite three examples: 

• One favoured case is the fi~+^'^C{0^) —>■ ^^6(1+) ± transition. In a ^ 1 
capture, there are three independent observableJ^: viz. the capture rate, T^'"^^, the 
polarization, Pav, of the recoil nucleus (here ^^B) along the muon spin direction, 
and the longitudinal polarization of the recoil nucleus, Pl, where 



P,, = {3-(T^)/J, 
Pl = (J-v)/J, 



(102) 



and (cr^) = is the muon polarization, J the nuclear spin, and v the recoil 
direction. Since the helicity of the muon neutrino, h^^, is — ^, ^^B can have only 
m = — lorm = spin directions in the recoil frame. The helicity of the final state 
(^^B+z/^j) is thus — I or for which the amplitudes are denoted by M_i/2 and 
M+i/2 respectively. The observables are expressed in terms of these amplitudes 
by 



.cap _ C{2 + x^ 

> 2, 



Pl = -2/(2 + (103) 

where x = a/2M_|_i/2/M_i/2, and C is a rate constant depending on known quan- 
tites. An interesting feature is that x can be determined from the ratio of Pav to 
Pl- 

i? = Pav/PL = -i(l + 2x), (104) 

when = 1. This ratio is not only more sensitive to x than either observable 
separately, but it is virtually immune to the systematic uncertainities associated 



with absolute measurements of recoil polarizations. Roesch et al. tl3 determined x 
in a novel recoil-implantation experiment that enabled Pav and P^ to be measured 
simultaneously. The parameter x depends not only on the pseudoscalar coupling 
constant, but also on nuclear matrix elements and other weak-interaction form 
factors, such as weak magnetism. With shell-model calculations of nuclear matrix 
elements, and with beta-decay experiments to establish the form factors (see 



sect. p.3[), the dimensionless quantity rrifj^g^/ is determined to be 9.0 ± 1.7. This 
result is slightly larger than the PCAC prediction, Eq. but is consistent with 
the RMC results in Table 7. 

Another interesting case is the /i~-|-^^O(0^) — > ^^N(O^) + z/^ transition and its in- 
verse, the beta decay ^^N(0~) "'^^O(0^)+e~+I7e, for which both transition rates 
have been measured The analysis here depends on the shell-model calcula- 



tions used to compute matrix elements of parity-changing operators. Gagliardi et 



alx^ obtained m^gp/gj^ = 11±2; however, subsequent evaluations by Warburton 



et al. — and by Haxton and JohnsonEH, with much improved shell-model wave 
functions, resulted in smaller values: in the first case, mp^gp/ g^ = 7.4 ± 0.5 and, 
in the second, it was in the range 7-9. Both are consistent with PCAC. 

Lastly, Gorringe et al^^ recently reported obtaining a value of m^gp/gp^ from 
measurements of the rates and hyperfine dependencies of yU~ capture on ^^Na to 
explicit daughter states in ^^Ne. The IS* atomic ground state of a muonic atom 



in the case of a Jj 7^ nucleus, is split into two hyperfine (HF) states with 
angular momenta = Jj + 1/2 and F_ = Ji — 1/2. For examples where the 
transition decreases the spin of the nucleus by one unit {e.g. |^ —>■ in fi~ ^^Na), 
capture from the F_ state is generally dominated by Qj^ and the Gamow- Teller 
(GT) matrix element, whereas capture from the F+ state is generally dominated 
by gp and the same GT matrix element. Consequently, the ratio of /^-capture 
rates, A+/A_, is a strong function of gp/gx and, to the extent that the GT 
matrix element dominates other terms, is independent of nuclear wave functions. 
Hyperfine transitions lead to the depopulation of the upper hyperfine state of 
the /i~ atom (F+ in ^^Na) in favour of the lower hyperfine state and yield a 
li~ capture time dependence that is sensitive to A+/A_. Gorringe et a/.t23 were 
able to measure this time dependence in the /i^ capture on ^^Na to the 1017 
keV and 1833 keV I"*" states in ^^Ne to obtain m^gp/gp, = T.Qlil'l, and < 7.1, 
respectively, consistent with the PCAC prediction. 



In summary, our knowledge of gp in nuclei is very sketchy and depends heavily 
on nuclear-structure evaluations for its determination. There is a hint in light nuclei, 
^^C and ^^O, that gp might be slightly enhanced over the PCAC value, while in heavy 
nuclei it seems almost completely quenched. Some support for the first remark comes 



from a calculation by Kolbe, Langanke and Vogelt^j, who with the continuum random 
phase approximation are able to reproduce the total capture rates for ordinary muon 
capture on the light nuclei, ^^C, ^^O and ^''Ca, without renormalization of gA or gp. 
The calculated rates would be significantly lower if in-medium quenching of coupling 
constants were introduced. 

That gp might be significantly renormalized in heavy nuclei was first discussed 
by Ericson et al. ^^^^^^^^^ who link the pseudoscalar coupling term to the pion-pole 
graph, described by Eq.|6l|. In the model of ref.^ the quenching is evaluated to be 
fi'p.off = 0.3(7p for infinite nuclear matter. This work is based on only p-wave pion-nucleon 
interactions. Recently, Delorme and EricsonE^l pointed out that the repulsive s-wave 
interaction significantiv alters this estimate, and their revised value is now gp^^ff = 0.6gp. 



Kirchbach and RiskaH3 also stress that the effective induced pseudoscalar coupling 
constant in nuclei originates from four different sources: the space-like and time-like 
parts of the axial current and the space-like and time-like parts of the pion pole, each 
of which is renormalized in a different way in the nuclear medium. The net effect for 
radiative muon capture in heavy nuclei in their calculation is a substantial quenching 
of transition rates, which increases with mass number as suggested by the experimental 
data. 

4.4- Meson- exchange currents in time-like axial- charge transitions 

In contrast to the quenching evident in allowed Gamow- Teller transitions dis- 
cussed in sect.[4.1|, transitions originating in the time part of the axial- vector current 



appear enhanced. Again there are two sources of explanations, both of which contribute 



to the final outcome. First, despite use of tlie best available shell-model calculations, 
the shell-model configuration space has inevitably been truncated, and corrections for 
this truncation are necessary. These are the core-polarization corrections discussed in 
sect.[4.1|. Second, the one-body operator from the impulse approximation is augmented 
by meson-exchange currents in which the axial-vector current prompts a meson to be 
emitted by one nucleon and reabsorbed on a second nucleon. Such processes lead to 
two-body operators. The pion, as the lightest-mass meson, produces a long-range op- 
erator and so it is the most important meson to consider for exchange currents. In any 
evaluation of these currents, the pion-absorption amplitude is well known, being deter- 
mined from the long-range part of the nucleon-nucleon force, but the pion-production 
amplitude is less certain. In beta decay, the momentum transfer is a few MeV and 
is certainly small relative to the nucleon rest mass. If, further, it is assumed that, in 
describing the pion-production amplitude, it is sufficiently accurate to take the limit 
that the pion momentum is small as well - i.e. the "soft-pion limit" - then this ver- 
tex is determined in a model-independent way through the exploitation of low-energy 
theoremsE£§. The net result is that the pion-production amplitude is given by the com- 
mutator, (1//^) [Qs, A^], where is the hadronic axial current, the axial-charge 
operator and the pion-decay constant. Were we dealing with zero-mass quarks with 
point couplings, then the evaluation of this commutator would be trivial: the charges 
would be generators of closed SU{2) x SU{2) algebra and the commutator would be 
determined from the group algebra. These same commut ators are nonetheless adopted 



in this context, following the premise of current algebrali2l't23 that the currents given 
by the quark model exhibit algebraic properties that are of more general validity. From 
the SU (2) X SU (2) group algebra, the commutator is 

[g^,4] = ze,,,\/;, (105) 

where the superscripts are isospin indices and Cabc the completely antisymmetric tensor 
in these indices. What this equation expresses is that the leading term in the weak 
pion-production amplitude is a current proportional to a vector current, V^. 

To continue, we evaluate this commutator between nucleon spinors and take the 
nonrelativistic limit. For the case in which the Lorentz index, fi, is time-like the com- 
mutator is of order 0{1), while when fi is space-like, it is of order 0{p/M). Here p is 
the momentum lost by the nucleon making the transition, and is small by comparison 
with M, the nucleon mass. By contrast, an evaluation of the axial current, A^, itself be- 
tween nucleon spinors - the impulse approximation - leads, in the nonrelativistic limit, 
to terms of order 0{p/M) when fi is time-like, and of order 0{1) when fi is space-like. 
Thus, from these limits, it is evident that meson-exchange currents (MEC) will be more 
important in axial-charge transitions (originating in the time component of the axial 
current) than in Gamow- Teller transitions (originating in the space component). This 



observation was first stressed by Kubodera, Delorme and RhoE£; . 



In sect. 4.1, we discussed corrections to some allowed Gamow- Teller matrix ele- 



ments and from Table 6 we see that the MEC correction is small: of the order of one 



percent. The isobar contribution is not part of the soft-pion MEC correction. This is 
because the pion-production amphtude at the vertex where a nucleon transforms to 
an isobar is of order 0{q/M), where q is the pion momentum, and vanishes in the 
soft-pion hmit. However, for the space part of the axial current, the soft-pion terms are 
suppressed - of order 0{p/M) - and are comparable in size to the correction terms. In 
this situation the soft-pion theorems by themselves are of little use. Rather, one must 
start with a model Lagrangian and explicitly evaluate the meson-exchange currents it 
would generate. 

For the rest of this section, we will only consider the time part of the axial 
current, for which the soft-pion terms are not suppressed: this leads to what are called 
axial-charge transitions. The two-body operator, after the nonrelativistic limit has been 
taken, is Fourier transformed to coordinate space and the rank-0 multipole is projected 
out, givingS'^nl 

Ao(soft-vr) = -^f^'^'^^-^YiM^i'Ti X ^2), (106) 

where cr+ = cr-^ + cr^, r = ri — r2, and Yi{x) = (1 + l/x)e~^ /x, with x = rriT^r. The 
corresponding one-body operator from the impulse approximation is 

Ao(IA) = Qa^^t. (107) 

The interesting question is how big are the matrix elements of the pion-exchange 
terms relative to those from the impulse approximation. This, then, requires a nuclear- 



structure calculation for the nucleus in question. Delormell^j has given an estimate for 
nuclear matter and finds the ratio y4o(soft-7r)/y4o(IA) to be large, being of order 60% at 
nuclear-matter densities. Axial-charge transitions are clearly heavily renormalized by 
pion-exchange currents. This prediction, which comes from soft-pion theorems and cur- 
rent algebra, is unfortunately hard to test experimentally because such transitions are 
difficult to isolate. There are two possibilities: one can either examine the small 1/M 
recoil corrections in an allowed transition, or one can look at transitions in forbidden 
decays. We will consider examples of both. 



4.4.1. Recoil corrections in allowed transitions 

The A=12 triad has long been the testing ground of nuclear beta-decay theories. It 
has been used to verify the conserved vector current (CVC) hypothesis through the 
identification of weak magnetism terms, and it has ruled out any possibility of size- 
able second-class currents (see sect.|2]^). By yielding the 1/M recoil correction in an 
allowed Gamow- Teller transition, it has also orovided a determination of the axial- 
charge matrix element. Various experiment J^@ii00Ei'Ei have measured the corre- 
lation between nuclear alignment and the direction of the emitted electron following 
nuclear beta decay from oriented ^^B and ^^N to the ground state of ^^C (see Table ^). 



The electron's angular distribution can be represented by a Legendre series, and the 
coefficient of the P2{cos9) Legendre polynomial is the quantity that must be measured. 
If we denote this coefficient as a+ for the positron emission from ^^N and as a_ for 
the electron emission from ^^B, the difference, a_ — a+, isolates the weak magnetism 
matrix element (together with the second-class current term, if present) while the sum, 
a_ + a+, isolates the axial-charge matrix element. It is the ant i- Hermit ian nature of the 
axial-charge operator that produces this sign change. The average of the experimental 
alignment data in Table ^yields + = —(2.823 ±0.078)GeV~^. Theoretically, the 
combination a_ + q;+ is related to nuclear-structure matrix elements via 



I 4 1 

a.- -\- — ~3 2m2/ 

y = l + 2{r{cT-V))/{(T), (108) 

where the matrix element (r(cr • V)) arises from a projection of the rank-1 multipole 
from the impulse-approximation axial-charge current. The experimental value for y is 
yexpt = 3.98 ±0.11.^ 



Morita et aliHH calculated these matrix elements in the impulse approxima- 
tion with shell- model wavefunctions that describe the ground states as Ohu {Op)~^ 
configurations in the closed-shell at ^^O. The results span y = 3.6 ± 0.1 depend- 
ing on the choice of residual interaction and are close to the experimental result. 
There seemingly is little evidence for a significant correction arising from two-body 



meson-exchange currents. Later, Koshigiri et alx^ realized that the shell-model space 



was quite restrictive and that enlarging the model space would significantly influence 



the (r(cr • V)) matrix element. Using perturbation theory, Koshigiri et a/.t23 demon- 
strated that this correction depended principally on the strength of the residual tensor 
force. Their final result is y = 3.6 — 1.2 -|- 1.3 = 3.7, with —1.2 being their core- 
polarization estimate with the Hamada- Johnston tensor force and +1.3 being the 
meson-exchange current enhancement. Guichon and SamourE^l, in a similar calcula- 
tion, obtained y = 3.0 — 0.4 + 1.0 = 3.6, where the smaller core-polarization corection 
of —0.4 arises with the Sussex interaction, which has a weaker tensor component. 
Clearly, details of the analysis depend on the choice of residual interactions, model 
spaces, wavefunctions, and so on, but if the predicted meson-exchange enhancement is 
present, it is in all likelihood masked by the core-polarization correction. 



4.4.2. First-forbidden beta decays 

In the case just covered, the axial-charge operator leads to a small 1/M correction 
in allowed decay, which has to be isolated from other small 1/M recoil corrections. 
By contrast, the time part of the axial- vector current contributes a leading term in 
first-forbidden weak transitions. Unfortunately, the space parts are not negligible: the 
familiar Gamow- Teller operator, retarded by one unit of orbital angular momentum, 



contributes with comparable magnitude. In first-forbidden decays, tlie impulse approxi- 
mation produces six one-body operators: one, the axial-charge operator cr • V/M, from 
the time part of the axial- vector current; three from the space part of the axial- vector 
current, [r, cr]^^\ with i? = 0, 1, or 2 being the tensorial rank; one from the space part 
of the weak vector current, V/M; and one from the time part, r. Angular momentum 
selection rules limit the number of these operators that can contribute. For example, 
in O"*" to 0~ transitions, only two rank-zero operators contribute; in 0^ to 1~ the three 
rank-one operators contribute; while in to ^ transitions, the five rank-zero and 
rank-one operators are involved. We limit our discussion to these three cases: thus, the 
rank-two operator will play no role. 

The notation for the first-forbidden beta-decay matrix elements is taken from 
Warburton et a/.B23'c21l and we define 



= g,{JfTf\ljjCT.VT\lJ,T,)C 

= -gAJfTf\\\a^rT\lJa^)C 

= g^V2 {Jf Tf III I a, rfW ||| J, T,) C 

Mf = -(J/ III rr III J, Ti)C 

Mf = -(J; III J^Vr III J, T,)C, (109) 

with C = {Jf /{\^Ji)){TiTzil =F 1 I TfTzf). The notation is the same as that intro- 
duced in Eq. ^ for allowed transitions. The conserved vector current hypothesis may 
be usedini to obtain an alternative expression for Mf in terms of M^: 

Mf = E^Mf , (110) 

where E^, is the photon energy of the analogous El gamma transition. In the Behrens- 
Biihringn23 formulation, which is being followed here, the beta-decay formulae are de- 
rived by expanding the lepton radial wavefunctions in powers of the lepton's mass and 
energy, and of the nuclear charge. In this treatment, inclusion in the radial integral of 
an extra factor that depends on the nuclear charge distribution introduces additional 
matrix elements obtained from Mq , and and denoted Mq', Mf and M"'. 

The ratios of the primed to unprimed quantities, denoted r'g, r'^ and r^, are insensi- 
lioo 



tivaiHU to details of the nuclear structure and are roughly 0.7. In the limit that small 
lepton-energy-dependent terms are dropped - the ^-approximation of Warburtonisil - 
the formalism for first-forbidden decays exactly matches that for allowed decays: viz. 



Eqs. 93 and 94. Thus 



6146 ± 6 



with fi(0) = I M(0) I 2, B{1) = I M(l) | 2, and 



M(0) = e^,,M^ + asM^ 

M(l) = auM^-a^M^, (112) 

where 



«5 = r'sC + lWo 

ax = E^-r'^^-lWo. (113) 

Here Z is the atomic number of the daughter nucleus in beta decay, Wq the end-point 
electron energy in electron rest-mass units, ^ = aZ/{2R) with a the fine-structure 
constant and R the radius of the nuclear charge distribution in electron Compton- 



wavelength units. The matrix elements given in Eq. |109| , and the a-factors in Eq. |113 
are computed in units that render them dimensionless. 

Lastly, we have included in Eg. |112| a factor, emec, representing the expected en- 



hancement of the axial-charge matrix element, Mq, over the impulse- approximation 
value due to meson-exchange currents. The procedure is to compute the matrix ele- 
ments with the best available shell-model wavefunctions, incorporate core-polarization 
corrections, and, by comparison of theory with experiment for several beta-decay tran- 
sitions determine the average value of the enhancement factor. This has been accom- 
plished by Warburton et al. for three mass regions. In the A = IQ region, the average^ 
from three beta decays and the /i~ capture on ^^O is emec = 1.61 ±0.03; in the A = 132 



region, the averagecSJ from three beta transitions is emec = 1.82 ±0.07; and for eighteen 
transitions in the A = 208 region, as re-analysed by Warburton and TowneiiH^, the 
average is emec = 1-79 ± 0.04. In the latter two cases, the results are dependent on the 
choice of the residual interaction used for the core-polarization correction. The values 
given are based on the Bonn nucleon-nucleon G-matrix, which has a weaker tensor com- 
ponent than, for example, the Paris nucleon-nucleon G-matrix. For the A = 208 region, 
the Paris G- matrix leads to a larger enhancement factor of emec = 1.91 ±0.04. These re- 
sults show that a substantial enhancement of the axial-charge impulse-approximation 
matrix element is required, and that there is some evidence for larger enhancement 
factors in heavier nuclei than in light nuclei. 



Shell-model calculations of emec have been given by TowneillHJ for closed-shell- 
plus-one configurations, for which the evaluation of the matrix element of the soft-pion 
operator of Eq. |106| involves the calculation of two-body matrix elements between va- 



lence and core nucleons and summed over all core nucleons. Such a calculation produces 
some mass dependence in emec because the number of core orbitals is changing from 
light to heavy nuclei, and the oscillator frequency parameter used in the radial wave- 
function decreases with increasing nuclear size. The results for the soft-pion operator 



are emec = 1-56 at A = 16, and emec = 1-76 at both A = 132 and A = 208. This is in 
reasonable accord with the values deduced from experiment. 



Calculations that go beyond soft pions, invoking heavy mesonsllHJ'c23 or chiral 



perturbation theoryc2ZI, and calculations that explain the mass dependence as 'Brown- 
Rho' scalingini can be examined in the cited references. 

5. Beyond the Standard Model 

5.1. Right-hand currents 

Maximal parity violation in weak interactions is built into the Standard Model 
as seems to be demanded by experimental evidence at low energy. At higher energies, 
one can envisage a restoration of parity symmetry through model extensions. In the 
charged-current sector, these models introduce additional, predominantly right-handed 
charged gauge bosons to complement the predominantly left-handed gauge bosons 
W^. The right-handed bosons acquire, via spontaneous symmetry breaking, a mass 
777.2 that is larger than the mass mi(~ 81 GeV) of the observed W"^. By construction, 
the left-handed couplings dominate at low energies with parity-symmetry restored at a 
higher-energy scale. The weak-interaction eigenstates - Wl, which corresponds to the 
Standard Model's W gauge boson, and Wr - are linear combinations of the physical 
bosons Wl and W2 with 



Wl = Wl cos C - Wr sin ( 

W2 = WLsmC + WncosC, (114) 
where C is a real mixing angle. The simplest of the model extensions is known as mani- 



fest left-right symmetryn22l, in which the physical left-handed and right-handed currents 
have identical coupling strengths and transformation properties, and one current is ob- 
tained from the other via 75 —75. In the pure V — A theory, the relevant part of a 
semi-leptonic weak interaction can be written in a schematic form as 

Tf, = ^Aiy- ^)] - ^)] = [VV + AA- {VA + AV)] , (115) 

where the {V — A) in the first bracket represents the hadronic current and in the second 
bracket the leptonic current. The modification to incorporate right-hand currents is 

T^, = £^^[ciV-A)-s{V + A)][c{V-A)-siV + A)] 



f-^ [s{V -A) + c{V + A)] [s{V -A) + c{V + A)] 

87772 



G 



^Vud [VV + r]AAAA + VAviVA + AV)] , (116) 



Table 8. Modifications to various observables of V — A theory through the introduction of 
right-hand currents evaluated to leading order'* in the small quantities C and 5. 



V - A Modification 



Fermi decay 


GIX'{M,)' 


G2(Mv)2(l + 2C-25) 
G2A2(M,)2(l + 6C-25) 
G2(l + 4C-25) 


GT decay 


Muon decay 


Fermi/muon decay 


G2(1 + 3A2) 


K?,(Mv)2(l-2C) 
Gl [(1 + 2C - 26) + 3A2(1 + 6C - 28)] 

o A2{1+(C+5)2)+A(1+C'+5^) 
^(l+(C-5)2)+3A2(l+(C+5)2) 


Neutron decay 


Neutron-decay asym., A 


9 A2+A 
^1+3A2 


Fermi long, polariz., P[ 


1 


1 - 2(C - 5)^ 


GT long, polariz., Pf^^ 


1 


l-2(C + 5)^ 


Ratio /PF 


1 


1 + 8C5 



"Note that radiative corrections are neglected, and Gv = G^Vud- 



where c is cos^ and s is sin^, and 

e^'o + 1 

= -^^TTT -1 - 2C + 25, (117) 
e^o + 1 

with e = (1 + tan(^)/(l — tan(^) and 5 = m\/rn^. In Table 8 we show the modifications 
to various observables in weak interactions evaluated to leading order in the small 
quantities C, and 5. Pure V — A theory is recovered in the limit that C and 5 — > 0. 

Beta asymmetry. As discussed in sect. ^]2|on neutron decay, the asymmetry in the 
distribution of electrons relative to a polarization in the parent neutron is governed by a 
parameter A that depends on the ratio. A, of axial-vector and vector coupling constants, 
A = G'^/G'y. With the inclusion of right-hand currents, A becomes a function of A, ( 
and 6 as given in Table 8. Likewise, the ratio of the neutron's ft value to that of the 
pure Fermi superallowed emitters is also given in terms of A, ( and 6. It is then a 
straightforward matter with these inputs to eliminate A and, allowing for experimental 
uncertainties, to set limit s^llS'lllill on ( and 6. A similar analysisiij can be performed 
for the mirror transition in the decay of ^^Ne, for which the beta asymmetry and lifetime 
have been determined. Carnoy et a/.S have also performed a combined fit to both data 
sets and obtain at the 90% confidence level: 5 = 0.12 ± 0.04 and C = -0.003 ± 0.009. 



The mass-ratio parameter 6 differs from the pure V — A result by nearly three standard 
deviations and this requires the mass of the predominantly right-hand boson to be in 
the range 200 to 300 GeV. 

Longitudinal polarization. Measurements of absolute electron polarizations in 
Gamow- Teller decays by van Klinkenill have a combined statistical precision of 0.9% 
and absolute accuracy of 1.6%. The most precise measurements are for highly hindered 
decays for which recoil and forbidden corrections might be important. However, if the 



evaluation of van KlinkenEla is taken at face value (Pf^ = 1.001 ± 0.012), a hmit of 
I ((^ + C) I ^ 0.10 is obtained at the 90% confidence level. 

An entirely different class of experiments is represented by two recent measure- 
mentsii^^ii^ of the ratio of positron polarizations / Pj^^ for pure Fermi and pure 
Gamow- Teller transitions. These are believed to be reliable measurements, with preci- 
sion of less than 0.5% for the polarization ratio. In addition to being free of absolute 
analyzing-power uncertainties in the polarimeter, the techniques employed in these 
measurements reduce many other systematic uncertainties. The results place a limit 
on the product of the two parameters: | C^l = 0.00015 ± 0.00034. 

CKM Unitarity. If the CKM matrix is taken to be unitary in three generations, 
then a value for V^d can be obtained from V^^ = 1 — V^^ — V^^^ with recommended values 
for Vus and V^t taken from ref.i. The result differs from the value determined from 
superallowed Fermi decays discussed in sect. ^.2.^1 (see also Fig. 2). If the difference 



is ascribed to right-hand currents, then a value is obtained for the mixing angle of 
( = 0.0018 ± 0.0007, which is two standard deviations from the pure V — A result. 

Muon Spectrum. Without radiative corrections, the muon differential decay rate, 
/x+ ^ e+ + z/g + z7^, integrated over e"*" spin directions, is given by 



d^r 



a;^dxd(cos i 



77? T — 1 

(3 - 2x) + (|p - l)(4x - 3) + 12--^— -^r? 



X 



(2x-l) + (|5^-l)(4a;-3) ^P^cosO. (118 



Here, x is the reduced electron energy Ee/E^ax, where -Emax = {rnl + m^)/(2m^) = 
52.83 MeV is the maximum energy, and rrie and are the particle masses. The 
angle between the positron momentum and the muon polarization vector, P^, in the 
IJ,~^ rest frame is vr — ^. The muon-decay Michel parameters are p, rj, ^ and 6m, the 
values of which in a pure V — A theory are: p = 3/4, ?7 = 0, ^ = 1 and 6m = 3/4. 
The experiment of Jodidio et a/.Eli was designed to measure the positron spectrum 
for X > 0.85 and cos 6' > 0.95, which is sensitive to the combination of parameters, 
^P^6m/ p- This combination equals 1 for pure V — A and 1 — 2(25^ + 26( + C^) with the 
inclusion of right-hand currents. The measured lower limit^ii of 

CP^.6m/p > 0.99682 (119) 
sets limits on the right-hand-current parameters at the 90% confidence level of 5 < 0.04 
or ms > 406 GeV and -0.056 < C < 0.040. 



It should be noted that the hmits on right-hand currents from muon decay and 
/?- asymmetry measurements disagree with each other at the 90% confidence level. This 
discrepancy, however, disappears£i3 at the 3a level. 

5.2. Scalar Interactions 

As has already been noted, the V — A structure of the weak interaction in the 
Standard Model has been put in a priori to obtain agreement with experiments. Al- 
though the model agrees with all experimental results at presently available energies, it 
remains an interesting quest to search for deviations from the V — A structure. In the 
last section, we discussed right-hand currents and in this we consider possible scalar 
and tensor couplings. A more general form of the weak-interaction Hamiltonian was 
first written down by Jackson, Treiman and Wyldisl as 



Tfi — —7=yud 
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+ \i)p(T^ui^n iJe-(^^lu{CT + C'rp'^^)^^^ . (120) 

A pseudoscalar contribution has not been included, since it vanishes in leading order 
in nuclear beta decay. Further, the induced terms, such as weak magnetism of Eq. |T0|, 
are not included here. Lastly, time-reversal invariance is assumed so that the coupling 
constants Ci and C[ are real. The Standard Model limit is recovered by the replace- 
ments: Cv = C'y = g^, Ca = C'j^ = (7a and Cs = C'g = Ct = = 0. Note that this 
parameterization also includes right-hand currents when the possibility that Ci 7^ C'^ 
is admitted. In the nonrelativistic hmit, ippipn (My) and ippa^u'ipn — ^ (^a) in nu- 
clear beta decay, so the presence of scalar interactions leads to modifications in Fermi 
transitions, as tensor interactions lead to modifications in Gamow- Teller transitions. 

In 1984, Boothroyd, Markey and Vogel^ set limits on the strength of scalar and 
tensor interactions in a least-squares adjustment to 93 pieces of data from nuclear beta 
decay that do not depend on nuclear matrix elements. There are seven parameters to 
be adjusted if right-hand currents are permitted - namely Ci/Cy and CJCy - but 
just three parameters if C[ = Ci. At the 95% confidence level (2cr) their results for a 
seven-parameter fit are 



Cs 










< 0.23 




< 0.190 


Cv 




Cv 




Cf 










< 0.09 




< 0.085 


Ca 




Ca 





Cs + Co 



Cv 
Ca 



< 0.065 



< 0.01, 



(121) 



and for a three-parameter fit 



Cv 



< 0.0033 



Cj' 

c^ 



< 0.0011. 



(122) 



These limits on scalar and tensor interactions come mainly from measurements on the 
beta-neutrino angular correlation coefficient, a, and the Fierz interference term, b. If 
one writes 



Kij — CiCj + C^Cj, 
then these observables are defined in terms of the coupling constants as 



(123) 



with 
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e = {M^)\Ktt + Kaa) + {M^?{Kss + K 



VVj 
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The upper sign is for electron emission and the lower sign for positron emission. In 
what follows, we comment on developments since the Boothroyd et al. surveylHHl . 

j3-v angular correlations. The classic e — v correlation data were obtained by 
observation of the lepton-recoil effect on the energy distribution of stable daughter 
nuclei from ^He, ^^Ne, ^^Ne and ^^Ar decays. Such experiments are very difficult be- 
cause the recoil energies are only a few hundred eV. High precision results are only 
available for ^He decayi^i and neutron decayS. Neither of these are pure Fermi decays 
and therefore are not sensitive to scalar couplings. However, a recent study of ^^Ar 
(0^ — > O"*") and ^^Ar — >• i"*") decays by Schardt and RiisagercHJ measured with 
high precision the shape of the beta-delayed proton peak following the superallowed 
decay to its isospin analogue state in the daughter nucleus. In this case, the lepton 
recoil is delivered to a proton-unbound state, and the subsequently emitted proton 
has a resulting energy spread that exceeds the energy imparted to the recoil daughter 
by a factor of ~ 60. Because its particle decays are isospin-forbidden and the natural 
width of the proton-unbound level consequently so small, the shape of the proton peak 
predominantly reflects the e — v — "p triple-correlation asymmetry parameter, which 
for isotropic distributions equals the e — v correlation coefficient. Adelberger's analy- 
of the Schardt-Riisager dataEH^ improves the limit on the scalar constant from 



sisP^^ 



the seven-parameter ffi by a factor of two: 



C^ 
Cv 



< 0.12, 



C^ 
Cv 



< 0.10. 



(126) 



It is evident that isospin-forbidden /3-delayed proton spectroscopy provides a powerful 
probe of scalar couplings. 



Fierz interference term. A measurement of the Fierz interference term in a pure 
Fermi transition would also lead to a limit on scalar couplings. In sect. p.2| , a discussion 
is given on the superallowed Fermi emitters and a statistical rate function, /, is defined 
in Eq. ^ The presence of scalar couplings adds an additional factor to the integrand in 
Eq. ^of (1 +'yb/W), where 7 is (1 — (aZ)^)^/^ with a the fine-structure constant and 
Z the atomic number of the daughter nucleus. Integration over the electron spectrum 
with the introduction of electromagnetic corrections, as in Eq. ^ leads to the modified 
expression 

The presence of scalar coupling in pure Fermi superallowed transitions would therefore 
show up as a dependence of the JFt values on (VT"^) rather than being constant. A fit 
of the data in Table 2 to 7(1^"^) gives a result of 6 = -0.005 ± 0.007 or | 6 1 < 0.017 at 
the 90% confidence level. In the limit of no right-hand currents, b is given by —2Cs/Cv, 
with the minus sign for a positron emitter, and hence | Cs/Cy \ < 0.008. This result is 



less restrictive than earlier fitstJ'EH3EH3tHZE22 of this type, mainly because the current 



set of jFt values shows more dispersion than the earlier data. 

Commentary. Any evidence for right-hand currents or scalar interactions men- 
tioned in the last two sections should be taken with extreme caution. Generally, such 
evidence depends on a number of assumptions, some of which are not always justified. 
In the past, when an experimental anomaly has been reported it has usually receded in 
the light of subsequent, improved experiments. As Deutsch and Quiniil put it: "The 
Standard Model is discouragingly successful" . 
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